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Abstract 


The  feasibility  of  using  the  earth's  crust  under 
the  oceans  as  a communications  channel  for  very  low  fre- 
quency waves  Was  studied.  The  structure  of  the  crust  and 
its  electrical  properties  were  used  to  evaluate  the  conduc- 
tivity-depth profiles  found  in  the  literature.  Using  the 
impedance,  lE  /K  , a new  phase  function  was  defined  that 
not  only  made  the  phase  changes  at  the  boundaries  of  the 
waveguide  explicit,  but  also  allowed  the  use  of  a computer 
for  calculation  of  the  propagation  and  attenuation  constants. 

The  system  noise  was  assumed  to  come  from  two  main 
sources*  atmospherics  and  thermal  effects.  Calculation 
showed  that  beneath  the  ocean  the  atmospheric  noise  could 
be  neglected.  The  thermal  noise  within  the  waveguide  was 
calculated  to  be  -162  dBm.  The  attenuation  for  the,  TM^q  ^ 

mode  was  calculated  to  be  -1.35  dB  /km,  thus  allowing  a 
transmission  range  of  only  about  163  km  for  1 kW  of  trans- 
mitted power.  The  conclusion  was  that) despite  the  aesthetic 

. ' - 

appeal  of  the  idea^he  crust  under  the  ocean  was],  feasible 
as  a long-range  communications  channel  at  VLF  for  the  (TM^q 
mode . 
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VERY  LOW  FREQUENCY  PROPAGATION  IN  THE  EARTH'S  CRUST 


BELOW  THE  OCEANS 

I . Introduction 


Background 

Communication  systems  are  the  key  links  between 
command  centers  and  control  units.  Without  the  continuity 
of  uninterrupted  communications,  unity  of  command  and 
control  of  resources  cannot  be  achieved.  With  the  advent  of 
the  nuclear  environment  in  modem  wartime  situations,  the 
potential  total  disruption  of  conventional  radio  communica- 
tions due  to  ionospheric  disturbances  has  created  the  need 
for  a more  secure  form  of  communications  system  (Ref  1:20). 
Direct  electromagnetic  transmission  through  the  crust  pro- 
vides the  possilility  of  a system  that  v/ould  not  only  be 
independent  of  ionospheric  conditions  but  also  have  trans- 
mitter and  receiver  hardened  against  attack  as  an  inherent 
feature  of  the  system. 

In  1961  the  idea  of  using  the  earth's  crust  as  a 
waveguide  for  electromagnetic  waves  was  first  introduced 

(Ref  2).  Numerous  studies,  attempting  to  determine  the 

\ 

feasibility  of  using  such  a hypothetical  waveguide  for 
communications  purposes  have  been  published  in  the  succeed- 
ing sixteen  years.  However,  they  have  two  main  aspects  in 
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commoni  1)  they  all  deal  almost  exclusively  with  the  crust 
beneath  the  continents;  2)  they  conclude  that  although 
such  a system  is  not  precluded  by  present  data,  neither 
can  the  feasibility  of  such  a system  be  clearly  shown. 

Because  present  war  scenarios  envision  operations 
on  at  least  two  continents,  it  has  become  important  to 
study  not  only  the  feasibility  of  lithospheric  communication 
on  the  North  American  continent,  but  also  between  North 
America  and  Europe,  under  the  ocean.  Two  years  after  the 
publication  of  the  original  concept  of  lithospheric  comm- 
unication, there  appeared  the  only  article  found  that  dealt 
exclusively  and  explicitly  with  the  oceanic  crust  (Ref  3) • 
This  was  more  along  the  lines  of  a concept  paper,  and 
presented  no  theoretical  background  to  support  much  of  the 
data  presented.  Unfortunately  this  article,  like  the 
others  on  the  continental  crust,  was  inconclusive  in  its 
findings.  The  only  other  article  to  treat  the  subject  of 
suboceanic  lithospheric  propagation  mentioned  it  only 
passing,  as  a contrasting  situation  to  the  continental 
case  (Ref  ktjk).  One  possible  explanation  for  the  lack 
of  work  on  this  specific  problem  may  be  that  until  recently 
there  has  been  no  data  available  on  the  oceanic  crust, 
and  even  now,  there  is  little  data  for  depths  greater  than 
one  or  two  kilometers.  For  whatever  reason,  the  problem 
has  been  neglected  over  the  last  sixteen  years. 
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Problem 


The  primary  purpose  of  this  study  was  to  determine 
whether  the  oceanic  crust  can  sustain  a VLF  wave  over  a 
sufficient  range  to  make  lithospheric  communications  feasi- 
ble. The  VLF  and  lower  LF  range  was  chosen  because  this 
frequency  is  low  enough  to  permit  lithospheric  propagation 
to  some  extent,  but  high  enough  to  avoid  noise  from  the 

natural  ELF  fields  in  the  earth.  The  crust  under  the  ocean 

0 

was  modelled  as  a waveguide  of  infinite  horizontal  dimensions 
and  finite  vertical  dimensions.  Vertically,  the  waveguide 
was  defined  by  the  seawater-saturated  sedimentary  strata 
above,  and  by  the  Mohorovicic  discontinuity,  below.  Lateral 
discontinuities  in  the  lithospheric  strata  were  disregarded, 
as  was  the  curvature  of  the  earth.  Although  the  frequency 
range  for  VLF  is  defined  to  be  3-30  kHz,  some  work  was  done 
in  the  low  end  of  the  LF  range  (30-3OO  kHz)  (see  Figure  1). 
Since  the  ocean  of  concern  was  the  Atlantic,  the  desired 


Figure  1.  Scale  of  frequency  categories  showing  range 
treated  in  this  study. 


range  was  taken  to  be  on  the  order  of  4500  km  (about  2800 
mi) . 

The  conductivity  of  the  crust  was  modelled  as  a con- 
tinuous but  rapidly  varying  fvinction  of  depth.  The  electro- 
magnetic waves  within  the  guide  were  assumed  to  be  planar 
TM  waves,  travelling  in  a direction  parallel  to  the  waveguide 
boundaries  (see  Figure  2).  The  waveguide  had  a depth  on  the 
order  of  one  wavelength.  The  lower  boundary  of  the  guide  was 
defined  electrically  by  a rise  in  conductivity  due  to 
thermal  effects  of  the  higher  temperatures  deeper  within  the 
crust.  Although  not  documented,  this  rise  in  conductivity 
seemed  often  to  coincide  with  the  geological  boundary,  the 
Mohorovicic  discontinuity. 

The  total  lithospheric  communications  system  problem 
can  be  divided  into  two  parts*  the  antenna  problem  and  the 


Figure  2.  Waveguide  model  of  the  crust  and  TM  wave. 
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propagation  problem.  A detailed  analysis  of  the  transmission 
and  receiving  antennas  would  be  futile  if  propagation  were 
not  feasible  through  the  mediiim  between  them.  Therefore 
this  study  treated  only  the  propagation  problem.  The  key 
factors  were  the  propagation  and  attenuation  constants 
derived  from  electromagnetic  waveguide  theory.  These  con- 
stants were  used  to  calculate  the  transmission  power 
necessary  to  make  reception  possible. 

The  results  of  the  attenuation  calculations  were,  by 
their  very  nature,  heavily  dependent  upon  the  peculiarities 
of  the  particular  conductivity  profile  used.  Because  of  the 
lack  of  dependable  conductivity  data  for  the  deeper  regions 
of  the  oceanic  crust,  the  waveguide  model  was  developed  as 
independently  of  the  specific  conductivity  profiles  as  was 
possible.  The  objective  was  to  enable  future  researchers  to 
simply  apply  newly-found  data  to  the  model  developed  here, 
so  that  it  would  not  be  necessary  for  the  waveguide  theory 
to  be  reworked  just  because  the  data  was  changed. 

Approach 

This  problem  involved  the  use  and  evaluation  of  data 
and  techniques  in  three  different  fields  of  study*  geophysics, 
electromagnetics,  and  communications.  Because  of  this 
interdisciplinary  nature,  this  report  will  discuss  these 
topics  in  more  detail,  and  often  on  a simpler  level  than  is 
usually  encountered  by  the  individuals  who  are  specialists 
in  the  various  fields.  The  intention  was  to  attempt  to  give 
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the  reader  a comprehensive  view  of  the  problem,  regardless 
of  the  profession  in  which  she  or  he  is  a specialist.  The 
report  also  attempts  to  enable  the  non-specialist  to  read 
only  the  less  technical  portions  and  still  obtain  a compre- 
hensive sximmary  of  the  research  performed. 

The  report  is  divided  into  four  main  sections: 
geophysical,  electromagnetic,  communications,  and  data.  The 
Geophysical  Analysis  chapter  discusses  the  various  aspects 
of  the  crust  that  affected  the  problem:  its  structural, 
thermal,  and  conductive  properties.  The  Waveguide  Analysis 
chapter  derives  the  mathematical  tools  necessary  to  calculate 
the  attenuation  constants  and  describes  the  computer  method 
used  to  perform  calculations.  The  Noise  Analysis  chapter 
presents  the  basic  theory  of  atmospheric  and  thermal  noise, 
and  their  effect  on  the  communications  system.  Finally  the 
Data  chapter  details  the  precise  numbers  used  in  calculations, 
and  the  results  obtained.  Entire  books  can  and  have  been 
written  about  any  one  of  these  topics j this  is  an  attempt  to 
present  the  problem  from  a comprehensive  perspective. 
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Introduction 

The  problem  of  the  feasibility  of  a lithospheric 
communication  system  was  studied  using  two  types  of  models. 

The  mathematical  model  of  the  electromagnetic  waveguide  is 
discussed  in  a separate  chapter.  The  geophysical  models  will 
be  treated  here.  The  importance  of  geophysics  to  this 
problem  should  not  be  overlooked.  Only  with  a thorough 
understanding  of  the  earth's  crust  can  the  basic  conductivity 
and  thermal  data  be  derived  and  applied  to  the  electromag- 
netic propagation  model.  Unfortunately,  a review  of  the 
literature  revealed  that  the  electrical  and  thermal  proper- 
ties of  the  crust  are  the  source  of  much  controversy  within 
the  geophysical  community.  For  this  reason,  several  hypotheses 
concerning  factors  affecting  the  data  obtained  from  actual 
experiments  will  be  discussed.  A list  of  geologic  terms  has 
been  provided  in  the  Appendix  for  the  convenience  of  the 
reader.  It  should  be  noted  that  the  theories  discussed  are 
highly  simplified  here,  in  the  interests  of  providing  back- 
ground in  the  geophysical  aspects  of  the  problem  without 
overwhelming  the  reader  with  technical  details. 

The  primary  motivation  of  this  study  was  the  assumption 
that  the  oceanic  crust  was  significantly  different  from  the 
continental  crust.  Several  feasibility  studies  have  been 
done  on  lithospheric  communication  under  the  continents. 
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Only  one  conceptual  paper  has  discussed  the  subsea  crust, 
published  in  19^3  (Ref  3) • The  continental  studies  were 
published  since  then.  By  examining  the  structures  of  the 
subcontinental  crust  and  the  subsea  crust,  this  chapter  will 
show  the  differences  between  them  and  give  the  background 
necessary  to  understand  the  conductivity  profiles  used. 

Crust  Structure 

One  theory  of  overall  earth  structure  and  crust  form- 
ation has  been  illustrated  in  Figure  3*  This  theory  considers 
the  structure  of  the  earth  in  terms  of  a very  slowly  circu- 
lating system.  Large  convection  cells  have  formed  between 
the  central  core  of  the  earth  and  the  outer  crust,  within  the 
mantle.  Circulation  is  the  product  of  both  thermal  convection 
and  gravitational  forces.  The  hotter,  less  dense  rock 
material  rises  to  the  surface  while  the  cooler,  denser 
material  is  pulled  down  toward  the  central  core.  The  crust 
is  the  cooled  shell  on  the  outside  of  the  earth,  which 
emerges  and  is  reabsorbed  in  the  circulation  process  (Ref  5*228). 

This  theory  of  internal  earth  processes  is  highly 
consistent  with  another  theory,  which  models  the  crustal 
processes.  The  theory  of  plate  tectonics,  known  popularly 
as  "continental  drift,"  is  based  on  the  assumption  that  the 
earth's  crust  is  conposed  of  a number  of  solid  plates,  which 
"float"  on  the  more  fluid  mantle.  The  plates  are  theorized 
to  move  quite  slowly — on  the  order  of  about  one  centimeter 
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per  year — from  the  eruption  zones  to  the  subduction  zones. 

The  former  corresponds  to  the  areas  where  two  convection 
cells  meet,  both  pushing  material  to  the  surface,  whereas  the 
latter  corresponds  to  areas  where  both  convection  cells  pull 
material  from  the  crust  down  to  the  mantle,  towards  the 
central  core  of  the  earth. 

The  eruption  and  subduction  zones  model  is  based  on 
evidence  found  in  some  of  the  oldest  rocks  found  in  the 
crust.  The  exact  explanation  is  still  quite  controversial; 
however  the  above  models  serve  to  aid  in  londerstanding  the 
nature  of  the  crust's  structure.  Eruption  zones  have  been 
observed  to  correspond  to  oceanic  ridges  and  volcanic  island 
arcs.  The  subduction  zones  often  correspond  to  continental 
rift  valleys  and  mountain  ranges.  Thus  in  the  oceanic 
areas,  the  crust  is  relatively  thin,  due  to  the  "pulling 
out"  action  of  the  spreading  seafloor.  In  the  continental 

areas,  where  crustal  material  is  pulled  down  into  the  mantle, 

/ 

/ 

a corresponding  "piling  up"  of  crust  material  occurs, 
producing  a much  thicker  crustal  structure  under  the  conti- 
nents (Ref  6 1 228, 289). 

Because  of  extensive  seismic  and  compression  wave 
studies,  there  is  general  agreement  on  the  composition  and 
structure  of  the  crust.  The  three  basic  rock  types  are 
sedimentary,  formed  from  materials  deposited  primarily  by 
water;  igneous,  formed  by  lava  flows  of  magma  from  the 
mantle;  and  metamorphic,  formed  from  sedimentary  or  igneous 
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rock  that  has  been  radically  changed  in  form  by  pressure, 
temperatures,  and  other  forces,  deep  in  the  earth.  The 
crust  can  be  divided  into  continental,  oceanic,  and  transi- 
tional areas.  The  transitional  areas  refer  to  the  continental 
margins  both  above  and  below  sea  level  (Ref  7:238). 

The  upper  layer  over  most  of  the  crust  is  composed  of 
sedimentary  strata  with  an  average  depth  of  1-2  km.  Under 
the  continents  these  strata  may  be  as  thick  as  10  km;  under 
oceans,  they  may  be  altogether  absent  in  places.  However  the 
usual  depth  under  oceans  is  about  1 km  (Ref  3*32^).  Below 
the  sediments  on  the  continents  is  located  a layer  of  granite 
that  becomes  increasingly  basaltic  w'th  depth.  In  some 
areas  the  branite  and  basalt  are  separated  by  a sharp  break 
call  Conrad's  discontinuity.  The  granite  layer  may  be  as 
thick  as  20  km,  while  the  basaltic  layer  extends  to  the  lower 
limit  of  the  crust,  defined  by  the  Mohorovicic  discontinuity. 
In  continental  regions,  the  Moho  is  at  an  average  depth  of 
35-^0  km,  but  vinder  large  mountain  ranges  it  may  lie  as  deep 
as  80  km  (Ref  8:8).  Undersea,  the  thin  sediments  rest  on  a 
layer  of  basic  rock,  which  often  has  little  granite  in  it, 
but  resembles  the  base  of  the  continental  crust  in  its 
basaltic  nature  (Ref  3:324).  The  Moho  under  the  ocean  is 
at  an  average  depth  of  10-15  km  (Refs  3«324,  5*199) • The 
transitional  areas  have  characteristics  of  both  oceanic  and 
continental  regions.  All  three  general  crust  structures  are 
depicted  in  Figure  4. 


Figure  4,  Characteristic  depth  structure  for 
continental,  transitional,  and 
oceanic  crust  regions  (Ref  3*324), 


Little  is  kicwn  about  the  Mohorovicic  discontinuity  or 
the  nature  of  the  rock  around  it.  The  Moho  is  characterized 
by  a significant  change  in  the  velocity  of  seismic  waves  that 
cross  it.  One  theory  is  that  the  rock  material  undergoes  a 
phase  change  at  this  depth,  due  to  pressure  and  temperature 
effects.  Another  is  that  there  may  be  an  actual  chemical 
change  in  the  rock  at  that  depth  (Ref  5s 198). 

The  waveguide  for  lithospheric  communication  is 
bounded  by  the  sedimentary  strata  above  and  the  Moho  below. 

In  the  continental  structure,  the  dielectric  center  is  then 
primarily  granitic  to  vitreous  basaltic,  whereas  in  the 
oceanic  structure  it  is  crystalline  basaltic  (Ref  3:324). 

If  the  characteristics  of  silicon  dioxide  (quartz)  are  assumed 
to  be  similar  to  granite,  of  which  quartz  is  a major  component, 
then  it  is  interesting  to  compare  the  dielectric  constants  of 
quartz  in  vitreous  and  crystalline  form.  The  relative 
dielectric  constant  of  fused  quartz  is  about  3-75  (Ref  9: 

F-55) » while  for  the  crystalline  form  it  is  about  4.5  (Ref 
10:9-111).  This  significant  variation  in  dielectric  proper- 
ties is  a good  reason  to  reject  the  extension  of  continental 
data  to  the  oceanic  case.  Additionally,  the  dielectric 
difference  between  granite  and  basalt  are  considerable. 

Where  granite  contains  coarse-grained  quartz,  feldspars,  and 
micas  (Ref  11:166),  basalt  contains  only  feldspars  and 
pyroxenes,  which  are  in  a finer-grained  form  (Ref  11:31). 

Thus  there  appears  to  be  a substantial  difference  between 
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the  proposed  waveguide  media  under  the  continents  and  the 
oceans . 

Most  of  the  information  on  the  crust  has  been  deduced 
from  indirect  measurements  using  seismic  (earthquake)  waves 
and  soundings  of  compressional  waves  produced  by  explosions. 
Thus,  although  there  is  little  direct  evidence  of  the 
precise  crustal  structure,  the  general  outline  given  above 
was  assumed  to  be  the  structure  of  the  waveguide  considered. 

The  upper  boundary  was  formed  by  the  seawater-saturated,  thin, 
sedimentary  layer,  while  the  lower  boundary  was  formed  by  the 
Moho  and  the  upper  mantle  composed  of  dunite.  The  basaltic 
region  between  acted  as  the  proposed  waveguide  medium. 

Thermal  Characteristics 

The  thermal  properties  of  the  crust  were  important  to 
this  study  for  three  main  reasons:  1)  one  of  the  driving 
forces  of  the  convection  cell  system  was  thermal  circulation. 
Warmer  substances  have  a tendency  to  expand  and  their  decreased 
density  gives  them  a buoyancy,  making  them  rise,  while  cooler 
substances  contract  to  a denser  state,  making  them  tend  to 
the  bottom  of  the  system.  2)  The  lower  boundary  of  the 
waveguide  was  electrically  defined  by  an  assumed  rise  in 
conductivity  in  the  crustal  rock.  This  rise  in  conductivity 
(discussed  in  the  section  on  conductivity  profiles)  was 
attributed  to  the  rise  in  electronic  conduction  due  to 
higher  temperatures.  3)  A component  of  the  noise  power  in 
the  overall  lithospheric  communication  system  is  due  to 
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thermal  effects  within  the  system  ( this  will  be  discussed  in 
the  chapter  on  noise).  The  temperature  of  the  crust  under 
the  ocean  has  been  estimated  with  the  aid  of  geothermal  heat 
flow  measurements  at  the  surface  of  the  crust,  along  with 
estimates  of  thermal  conductivity,  radioactivity  in  the 
crust  and  mantle,  and  frictional  heating  in  the  crust  (Ref 

I 12:227). 

i 

i Heat  Flow.  A typical  temperature-depth  curve  for  the 

; earth  is  shown  in  Figure  5*  The  concave  shape  of  the  curve 

; is  due  to  the  addition  of  three  primary  heat  sources  in  the 

earth:  radioactive  elements,  friction,  and  pressure.  The 
principal  source  of  heat  is  the  decay  of  naturally  radio- 
active elements  in  the  upper  mantle  (ref  13s211).  Another 
|!  significant  cource  in  the  lower  crust  is  frictional  heat 

, produced  by  the  movement  of  tectonic  plates  across  the 

I 

mantle.  These  two  sources  add  t<k  produce  a very  steep 
temperature  gradient  in  the  first  200  km  of  depth.  Below 
this,  though,  the  much  less  steep  and  linear  increase  in 
temperature  due  to  pressure  effects  becomes  dominant,  since 
convection  and  radiation  are  minimal  (Refs  13«223,  14:264-6). 
Because  of  the  slightly  deeper  emplacement  of  the  heat 
sources,  due  to  the  thicker  overburden  on  the  continents, 
the  oceanic  temperature  gradient  is  less  steep  in  the  first 
100  km  than  the  continental  gradient  (Ref  13:220).  The 
continental  areas  of  crust,  in  addition  to  being  thicker 
than  the  oceanic  are  cooler,  due  to  their  being  older  rock. 

15 


Figure  5*  Thermal  profiles  for  typical  continental 

and  oceanic  heat  flows.  Dashed  line  shows 
assumed  melting  curve,.  (Ref  13*217) 

It  has  been  shown  that  seafloor  sediments  become  cooler  as 
they  age  (Ref  12i227),  and  this  appears  to  apply  to  conti- 
nents as  well. 

I Typical  heat  flow  rates  range  between  0.6  and  2.4 

2 

cal/cm  -s.  High  heat  flow  levels  at  the  surface  indicate 
an  area  of  geotechnic  activity,  since  these  are  also  areas 
where  the  hotter  mantle  material  approaches  the  surface. 

Friction.  In  addition  to  conduction  of  radioactive 
heat  in  areas  of  geotechnic  activity,  it  is  thought  that 
heat  is  produced  as  part  of  the  movement  and  reabsorption 
of  crustal  material  across  and  into  the  mantle.  It  has  been 
calculated  that  the  conversion  of  potential  energy  to  heat  in 
the  process  of  a plate  being  reabsopbed  into  the  upper  mantle 
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produces  a finite,  although  indefinite,  amount  of  heat  (Ref 
12i229}>  It  has  also  been  theorized  that  hot  areas  are  left 
in  the  wake  of  a plate  being  dragged  across  the  upper  mantle 
as  it  is  absorbed  in  a mountainous  region  nearby.  The 
resulting  frictional  heat  meltt-;  the  upper  mantle  and  produces 
a volcanic  island  arc  (Ref  12:228-30). 

Temperature  Profiles.  The  heat  flow  measured  ax  the 
surface  can  be  u.ed  to  find  the  temperature  gradient  by 
dividing  heat  floff  by  the  thermal  conductivity  of  the  crust 
material  (Ref  15*24).  However,  direct  measurements  of  thermal 
conductivity  extend  do’..-n  only  to  a depth  of  several  kilometers. 
Thus  belov;  this  depth,  thermal  conductivity  figures  are 
estimates,  based  on  other  indirect  measurement  methods.  This 
uncertainty  has  produced  some  disagreement  in  the  literature 
about  the  value  of  the  thermal  gradient.  As  shown  in  Figure 
5,  the  gradient  above  50  km  is  essentially  linear.  The  value 
of  the  gradient  given  here  was  15°/km.  Values  as  high  as 
20°/km  have  been  reported  for  the  oceanic  crust  (Ref  2:190). 
Continental  crust  gradients  are , for  the  reason  mentioned 
earlier,  lower,  although  they  also  are  as  high  as  oceanic 
gradients  in  some  areas.  Continental  thermal  gradients 
range  from  9°/kra  to  24°/km  (Ref  16*177). 

Conductivity  Profiles 

Factors  Affecting  Conductivity.  The  basis  for  the 
waveguide  attenuation  calculation  was  the  conductivity  profile 
of  the  oceanic  crust.  Because  the  conductivity  is  so 
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dependent  on  the  various  structure  factors  and  properties, 
it  was  necessary  to  examine  the  geologic  and  thermal  profiles 
before  studying  the  conductivity.  In  fact,  the  main  reason 
why  this  problem  is  still  being  debated  sixteen  years  after 
its  conceptualization  is  that  the  conductivity  profiles  cannot 
be  measured  directly  (or,  rather,  for  financial  reasons,  have 
not  been) , thus  leading  to  a disagreement  over  the  exact 
shape  of  them.  A review  of  the  main  factors  affecting 
conductivity  in  the  crust  will  make  the  uncertainties  more 
clear. 

There  are  two  basic  forras  of  conductiom  electro- 
lytic (ionic)  and  electronic.  The  former  refers  to  the 
flow  of  ionic  charge  carriers  in  a solution  of  electrolytes 
(dissociated  ions)  (Ref  17*11).  This  type  of  conductivity 
is  highly  dependent  upon  frequency  (Ref  18*348).  The  latter, 
electronic  conduction,  refers  to  the  flow  of  free  electrons 
either  in  metallic  ores,  or  in  other  ionic  substances,  when 
high  temperatures  have  already  broken  up  the  ions  Ref  17:11). 
The  basic  theory  of  lithospheric  conductivity  was  based  on 
the  combination  of  electrolytic  and  electronic  effects. 

The  conductivity  profile  is  affected  by  temperature, 
pressure,  and  water  content  of  the  rock  in  the  crust.  At 
temperatures  around  450°C,  it  is  theorized  that  electronic 
conductivity  effects  become  dominant  and  the  conductivity 
rises  at  a significant  rate  (Ref  4:8-9) . The  significance 
of  pressure  effects  on  the  conductivity  is  not  quite  as 
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well-defined,  since  pressure  effects  seem  to  counteract 
temperature  effects.  Assuming  a constant  density  of  crust 
mater^Al,  the  pressure 3 at  various  depths  may  be  calculated. 
Adding  to  the  density  pressure  is  pressure  from  stresses 
within  the  crust:  at  10  km,  the  density  pressure  is  O.27 
gigapascals,  while  the  combination  of  stress  and  density 
pressure  is  O.67  IPaj  at  30  km,  the  pressures  are  0.81  GPa 
and  1.2  GPa  (Ref  5*15)*  There  does  seem  to  be  evidence 
that  the  counteraction  of  pressure  effects  on  temperature 
effects  is  nontrivial  (Ref  19:158),  but  the  nature  of  the 
interaction  of  the  two  effects  has  not  yet  been  clearly 
def ined . 

The  other  major  factor  affecting  the  conductivity  is 
the  water  content  of  the  rock.  Water  occurs  in  three  main 
forms  in  mineral  structures:  1)  as  hydroxyl  ions,  locked 
into  the  molecular  structure  of  the  minerals?  2)  as  a free 
solution  in  the  pores  and  fractures  of  the  rocks?  3)  as 
films  in  capillary-sized  raicrocracks.  In  the  first  case, 
the  presence  of  absence  of  hydroxyl  ions  has  no  effect  on 
conductivity.  In  the  second  case,  the  water  serves  as  a 
weak  electrolyte  of  mineral  ions,  and  thus  raises  the  con- 
ductivity by  several  orders  of  magnitude.  The  effect  of 
interstitial  films  is  very  controversial?  however  there  are 
indications  that  they  also  raise  the  conductivity  of  the 
rock  by  several  orders  of  magnitude  (Ref  19:159-160). 

Methods  of  Measuring  Conductivity.  To  depths  of 
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about  3 Ion  I direct  measurements  have  been  made  on  the  litho- 
sphere, These  measurements  usually  are  made  in  deep  bore- 
holes produced  by  petroleum  exploration  wells.  Unfortunately, 
few  have  been  drilled  in  areas  with  thin  sedimentary  layers, 
down  to  the  granitic  and  basaltic  layers,  since  petroleum 
is  not  usually  found  in  these  areas.  Fewer  still  such 
boreholes  have  been  drilled  in  oceanic  areas.  There  are  also 
many  questions  about  interpretation  of  the  data  from  these 
boreholes,  due  to  such  factors  as  the  conductivity  of  the 
fluid  in  the  borehole  and  the  anisotropy  of  the  rock  forma- 
tions containing  the  borehole  (Ref  19*15^-155)* 

In  order  to  determine  conductivities  of  deeper  levels 
in  the  crust,  indirect  methods  have  been  developed.  There 
are  two  types*  measurements  on  the  crust  itself  and 
laboratory  measurements  on  samples.  The  former  category 
includes  soundings  using  magneto-telluric , four-electrode, 
long  antenna,  or  magnetic  loop  methods.  Magneto-telluric 
methods  involve  surface  measurements  of  the  natural  electro- 
magnetic field  along  a line  of  constant  magnetic  latitude 
(distance  from  the  north  magnetic  pole).  By  measuring  the 
field  fluctuations  and  subtracting  out  the  components  due  to 
external  cuases,  the  internal  fields  are  left,  giving  an 
indication  of  the  distribution  of  lithospheric  conductivity. 
Although  this  method  is  frequently  used  in  geophysical  pros- 
pecting, its  usefulness  in  conductivity  studies  is  not  clear, 
as  there  is  disagreement  as  to  how  the  data  should  be 
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interpreted.  In  any  case,  magneto-telluric  methods  provide 
maximum  values  for  the  conductivity  (Ref  20j405-407).  ’^‘he 

four-electrode  method,  one  of  the  oldest,  has  similar  problems 
to  the  magneto-telluric.  Long  antennas  and  magnetic  loops 
are  also  useful,  but  their  precesion  is  questionable.  All 
of  these  methods  are  useful  when  their  results  are  used  in 
conjunction  with  previous  data  on  the  geology  of  the  area 
and  other  techniques  of  measurement.  However  they  are  very 
difficult  to  use  alone . 

The  measurements  that  can  be  performed  in  the  labora- 
tory are  virtually  unlimited.  Samples  of  rock  from  drilled 
cores  or  rock  of  the  same  minerals  and  structure  as  those 
within  the  crust  can  be  tested  for  conductivity  under  a 
variety  of  pressure,  temperature,  and  humidity  conditions. 
However,  questions  have  been  raised  as  to  how  similar  these 
simulated  conditions  are  to  actual  crustal  conditions. 

One  of  the  reasons  that  doubt  has  been  cast  upon  all 
of  these  methods,  is  that  the  results  obtained  are  widely 
variant,  resulting  in  the  basic  controversy  of  the  litho- 
speric  communication  problem.  The  ideal,  of  course,  is  that 
assiaming  the  conductivity  profile  is  stable  over  a period 
of  time,  various  measurement  techniques  should  produce  conver- 
gent results  that  would  indicate  the  "correct"  values.  How- 
ever, the  obtained  values  diverge  substantially.  The  reason 
for  this  lies  in  the  nature  of  the  material  being  studied. 
Hocks  within  the  crust  are  composed  of  an  extremely  complex 
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three-phase  system:  the  mineral  matrix  which  is  porous,  the  ' 

( 

interstitual  solution  of  mineral  ions  and  water,  and  gas.  j 

i 

Thus  a laboratory  reproduction  if  these  exact  conditions  ' 

i 

is  difficult  if  not  nearly  impossible  to  achieve.  It  becomes  j 

even  more  difficult  to  interpret  measurements  taken  when  only  | 

some  of  the  lithospheric  conditions  are  reproduced.  Whereas  ] 

I 

measurements  taken  using  electromagnetic  soundings  are  not  i 

nearly  as  accurate  aind  do  not  have  nearly  the  resolving 
power  needed  to  determine  conductivities  accurately.  More- 
over, soiinding  methods  are  based  on  an  extremely  simple 
model  of  the  earth's  natural  electromagnetic  field--one  that 
ignores  cracks  and  fissures,  stratifications,  and  inhomogene- 
ities. Once  again,  it  is  very  difficult  to  interpret  the 
data  obtained,  even  with  the  knowledge  of  the  simplifications 
used  (Ref  3:20) . 

This  discrepancy  between  data  obtained  from  different 
sources  is  the  crux  of  the  controversy  over  the  feasibility 
of  lithospheric  commimication.  Assessment  of  this  feasibility, 
particularly  in  oceanic  regions,  was  made  difficult  by  the 
lack  of  data  and  the  conflicts  between  the  existing  data  on 
conductivities.  It  cannot  be  overemphasized  that  without 
reliable  conductivity  data,  any  calculation  of ' the  propagation 
of  electromagnetic  waves  in  the  crust  becomes  suspect.  The 
development  of  this  problem  seems  to  have  reached  the  point 
at  which  the  theory  has  been  developed  as  far  as  possible 
with  the  present  data.  For  further  advance,  new  data  will 
be  necessary. 
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Actual  Conductivity  Profiles.  The  lack  of  appropriate 
and  reliable  conductivity  data  became  a problem  in  this 
investigation  when  the  need  arose  for  a set  of  conductivity- 
depth  profiles  to  use  in  calculations.  There  were  only  two 
profiles  found  that  were  based  on  data  from  the  oceanic 
lithosphere  (see  Figure  6),  and  these  were  composites, 
derived  as  two  alternatives  from  the  same  data  base  by  the 
same  authors.  An  additional  drawback  of  these  profiles  was 
that  they  were  over  fifteen  years  old,  from  the  only  available 
attempt  to  discuss  the  specific  problem  of  oceanic  lithospheric 
communication  (Ref  3).  Despite  the  disadvantages,  these 
profiles  were  used  for  the  very  reason  that  they  were  the 
only  oceanic  profiles  found. 

The  other  five  profiles  were  chosen  for  comparison 
purposes.  As  profiles  of  continental  conductivity,  the 
attenuation  parameters  for  these  profiles  had  been  calculated, 
using  boundary  conditions  and  methods  appropriate  to  the 
continental  lithospheric  waveguide.  By  using  the  continental 
profiles  with  oceanic  boundary  conditions  it  was  hoped  that 
an  evaluation  of  the  effects  of  the  changes  on  attenuation 
due  solely  to  the  calculation  method  and  boundary  conditions 
could  be  made.  All  these  profiles  are  included  in  the 
Appendices. 

In  spite  of  the  widespread  controversy  about  the 
exact  mechanisms  and  values  of  the  conductivity-depth  profile, 
a general  profile  may  be  described.  Below  the  continents. 
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CONDUCTIVITY  IN  MHO/N 


the  conductivity  decreases  from  a relatively  high  value  due 
to  a high  water  content  in  the  rock  pores  and  cracks  (Ref 
4«8).  The  exact  reason  for  the  decrease,  down  to  a depth  of 
about  ten  kilometers,  is  a major  source  of  contention  in  the 
field  of  geophysics.  One  view  is  that  pressuie  closes  the 
cracks,  drying  the  rock  out,  and  thus  lowering  the  conduc- 
tivity. The  opposing  view  is  that  temperature  effects  are 
the  major  factor  in  decreasing  conductivity:  at  depths  where 
electrolytic  conduction  effects  start  being  curbed  by  electronic 
conduction  effects,  the  net  conductivity  decreases  (Ref  21:33^)* 
In  any  case,  it  is  agreed  that  the  conductivity  does  decrease 
down  to  a depth  at  which  the  electronic  conduction  dominates 
and  the  conductivity  rises  again,  due  to  rising  temperature. 

The  conductivity  of  the  subsea  crust  is  limited  at  the 
upper  surface  by  the  conductivity  of  sea  water.  The  water- 
saturated  sediments  have  a significantly  lower  conductivity 
which  then  holds  constant  or  decreases  slightly  with  depth, 
depending  upon  which  view  of  water  content  of  the  rock  is 
taken.  Then,  at  the  depth  where  the  temperature  increases 
to  a point  to  effect  electronic  conductivity,  the  conductivity 
increases  (Ref  12:230). 

In  order  to  understand  much  of  the  background  and 
controversy  concerning  conductivity  profiles,  it  was 
necessary  to  review  the  geologic  structure  and  geotechnical 
processes  in  the  crust.  The  calculation  of  the  attenuation 
of  the  electromagnetic  waves  in  the  crust  required  profiles 


of  the  conductivity  over  the  depth  of  the  crust.  Despite  the 
controversy  over  these  profiles,  the  ones  foimd  in  the 
literature  were  used.  Only  two  pertained  to  the  oceanic 
crust.  The  other  five,  taken  from  the  continental  crust, 
were  used  for  comparison  with  the  work  already  done  by 
other  researchers. 
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III.  Waveguide  Analysis 


Introduction 

Assumption  of  the  use  of  a vertical  monopole  trans- 
mitting antenna  led  to  the  assumption  that  only  TM  modes 
would  exi st  in  the  waveguide . The  waveguide  medium  was 
treated  as  flat  and  horizontally  homogeneous,  neglecting 
stratifications  and  discontinuities  as  well  as  the  earth's 
curvature . 

The  conductivity  within  the  crust  was  assumed  to  be 
significantly  less  than  the  boiandary  conductivities,  thus 
producing  a continuous  but  rapidly  changing  conductivity- 
depth  function.  For  this  reason  a phase  integration  technique 
was  chosen  over  both  boundary-value  matching  methods  and  WKB 
approximation  methods.  The  former  is  best-suited  for  a dis- 
continuous conductivity-depth  function,  while  the  latter 
applies  best  to  a slowly-varying  continuous  function.  Thus 
an  application  of  the  phase-integration  method  to  a simple 
horizontally  stratified  crust  v.'ould  require  some  modifica- 
tion to  account  for  the  extremely  rapid  changes  in  conductivity 
between  strata. 

Derivation  of  the  Seneral  Impedance  Equation 

The  electric  and  magnetic  fields  were  defined  as 

E(x,z,t)  = idcx-*^)  (1) 
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(2) 


ff(y,t)  = CeyHy(z)}  exp  i(  kx  - <ut) 

(Symbols  are  defined  in  the  List  of  Symbols,  p.  vii)  The 
exp  (-iwt)  time  dependence  will  be  suppressed  hereafter. 
The  conductivity  and  therefore  the  complex  electric  permit- 
tivity was  taken  to  be  depth-dependent; 


0-  = er{z) 

(3) 

€ = €( z) 

(4) 

= € • + -^ 

(5) 

Starting  with  Maxwell's  equations 

V X E - iwfilT  = 0 

(6) 

V X IT  + iu;e'E  = E 

(7) 

7"*  K IT  = 0 

(8) 

< 

• 

II 

o 

(9) 

and 

applying  the  vector  identity 

7 X (fA)  = 7f  X A + f(7  X A) 

(10) 

the 

following  relation  was  obtained; 

7 X E = 7x 

(11) 

- 

dE 

* * ike^E^(z)exp  i^  } x e„ 
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where 


7 X "ix  = V X e^  = 0 (13) 

and 

<l>  = kx  ( 14) 

Thus 

- - dE„ 

V X E = ey£  exp  iif  - ikE^exp  } (15) 

Substituting  into  Eq  (6) 
dE^ 

^ = ikE,  + (16) 

dz  z y 

To  obtain  a similar  expression  for  the  magnetic  field,  the 
identity  Eq  (lO)  was  applied  to  the  definition  Eq  (2),  and 
the  following  was  obtained: 

V X H = 7(H  exp  i4>)  x .e  (1?) 

dH 

= exp  i<>  + ike^H^exp  i<>  (18) 

Substitution  of  Eq  (18)  into  the  Maxwell  equation  (7),  when 
the  definition  for  the  complex  permittivity,  Eq  (5)i  was 
applied,  yielded 

_ dH 

7xff  + iw€E  = -ex-j“  exp  + ike^H  exp 

(19) 

+ i«^CEj^ex  + E„e^}  exp 

= 0 (20) 
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Thus,  Eqs  (19)  and  (20)  produced 

- _ _ _ 

^ * i“'CVx  ■"  ®i®z5  ' ° 


(21) 


So  that 


/ 


dH 


(22) 


and 


iffly  + iw^E^  = 0 


(23) 


since  the  coefficients  of  orthogonal  vectors  must  indivi- 
dually equal  zero,  if  their  sum  is  to  equal  zero.  Eq  (23) 
produced 


\ ' -1^  Hy 


(24) 


Substitution  into  Eq  (16)  yielded 

ifx  = icjtx  Cl  - _k^  1 

d z w2^g  ^ 

If  Eq  (22)  was  expressed  in  the  form 


^ = iu)€E, 

dz 


(25) 


(26) 


then  Eqs  (25)  and  (26)  were  the  two  basic  equations  to  be 
solved  for  the  fields  and  propagation  constant.  For  ease 
of  calculation  the  function  A was  defined 

E„ 


A = 


H. 


(27) 


Thus  the  derivative  became 
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I 


L 


^ = dA  u 4.  A y 
dz  dz  y ^ dz 


= iw  fiCl  - 


-/-3  H 


(29) 


where  the  definition  of  the  field,  Eq  (2)  was  applied. 
Substituting  Eq  (26)  into  Eq  (28),  the  equation  of  Eq  (28) 
to  Eq  ( 29 ) gave 


+ A(iu,€E  ) = icuCl  - -2^} 
Applying  the  definition  of  A,  Eq  (27),  the  equation 


( 30) 


«y  = iHiC  1 - «v 

^ 0)  Jlfc  J 


(31) 


was  obtained . Thus 


hI  = i‘-C-€A"  -H  fx  - 


(32) 


Normally,  this  equation  would  have  been  solved  by  direct 
integration.  However,  attempts  to  do  so  were  hindered  by 
the  presence  of  singularities  in  the  function  A( z) . These 
singularities  were  removed  by  setting 


A = -^  tan  u 

Substitution  into  Eq  (32)  produced 


(33) 


dA  _ . r 

1 2 du 

= T u ^ 


1^2 

- € tan  u + (i} 


(34) 


(35) 


where  the  latter  was  obtained  from  the  definition,  Eq  (33). 
Thus  the  final  equation  for  the  phase  function,  u,  was 


sin^  u - il  cos^  u { - wfx}  (36) 

This  was  the  eq.uation  sqlved  by  numerical  integration  tech- 
niques, 

Ntmierical  Method 

The  method  used  to  integrate  equation  (36)  and  obtain 
values  for  the  complex  propagation  constant  k had  two  primary 
functional  units*  a minimizer  and  an  integrator.  A detailed 
description  of  all  the  various  program  units  used  has  been 
placed  in  the  Appendix.  The  initial  and  final  boundary 
conditions  for  the  waveguide,  derived  in  the  following 
section,  were  input  to  the  minimizer,  along  with  an  arbi- 
trary initial  value  and  range  values  for  k.  The  integrator 
used  an  ei^th-order  predictor-corrector  method  to  calculate 
the  value  of  the  phase  function  at  the  lower  boundary  of  the 
waveguide.  This  value  was  compared  with  the  boundary  value 
given  to  the  minimizer.  The  minimizer  then  produced  a new 
value  of  k,  attempting  to  make  the  difference  between  the 
integration  results  and  the  boundary  condition  as  small  as 
possible.  The  successive  values  of  k were  based  on  random, 
gradient,  average,  and  jump  steps.  After  the  desired  number 
of  trials,  the  best  value  was  output,  giving  the  complex 
propagation  constant.  Both  the  minimizer  and  the  integrator 
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were  library  programs,  run  on  the  ASD  computer  system,  a 
CDC  CYBER  7^. 


Verification  and  Boundary  Conditions 

In  order  to  test  whether  the  computer  program  was 
actually  doing  the  desired  calculations,  three  special 
cases  of  the  waveguide  problem  were  used.  In  each  of  the 
three  cases,  the  bovindary  conditions  had  to  be  derived,  and 
the  values  of  the  propagation  constant  calculated  from 
known  waveguide  theory.  The  boundary  value  expressions 
were  used  in  the  computer  program,  and  the  resulting  k 
was  compared  with  the  theoretical  value.  In  all  three 
cases  the  agreement  was  within  2%. 

Case  I I Infinitely  Conducting  Boundaries  With  Perfect 
Dielectric . We  first  considered  the  case  of  a lossless, 
dielectric-filled  waveguide  with  two  parallel  and  infinitely 
conducting  boundaries.  As 


» 1 


b> 


(37) 


Substitution  in  Eqs  (25)  and  (26)  yielded 


dZ  |o_aO 


io.  E^}  = 


- ffE. 


(38) 


Sc 

dz  I 


If -i^oe 


1 w p <r 


= i«4CH 


y 
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(39) 


Eqs  (38)  and  (39)  were  expressed  in  matrix  form  and  solved 
as  an  eigenvalue  equationi 


r 


The  characteristic  equation  was  solved  for  the  eigenvalue  \ : 

2 

X+iu)H<r=0  (^1) 

\ = + ( -1  + i)  = ± (42) 


where  6 was  defined  as  the  skin  depth.  This  resulted  in  the 
eigenvector 


in  the  (-z)  direction.  Thus 

Re  ( X ) > 0,  X = ^ (46) 

was  chosen.  Thus  the  initial  condition  was 

= ^i^(i  - 1)  (47) 

When  the  conductivity  was  assumed  to  be  infinite  the  initial 
condition  tended  to  zero^  Using  this  as  an  initial  condition 
the  integration  of  Eq  (32)  was  performed.  Equivalently,  the 
initial  condition  for  u in  Eq  (36)  was  set  to  zero.  Since 
the  waveguide  medium  was  a perfect  dielectric,  the  dielectric 
constant  was  taken  as 

€ = €•  (48) 

since  the  conductivity  was 

<r  = 0 (49) 

However,  to  check  the  final  value  of  the  impedance 
function  at  the  lower  boundary  of  the  waveguide,  another 
derivational  method  was  used.  The  depth  dependence  of  the 
electromagnetic  fields  was  considered  as  a function  of 
exp  (+az),  where  it  was  desirable  to  obtain,  first,  an 
expression  for  a . The  magnetic  field  was  expressed  as 

IT  = (H^exp  i * H_exp  i ♦_)  9^  (50) 


where 
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kx  + a z 


(51) 


= 


From  the  identity,  Eq  (10)',  v x H was  calculated  using 


V(  exp  i <|)^)  = (exp  i iV 


(52) 


and 


V ♦ + = e k + e ar 


(53) 


Thus 


V X ± ©2^*  ^ * SyH^exp  i 

• i(e^k  + ©x®  ^ '*’  + 

Using  the  complex  form  of  the  dielectric  constant,  the 
Maxwell  equation  (7)  "became 


(54) 

(55) 


“ iwe 


V X H, 


Substitution  of  Eq  (55)  into  (56)  yielded 
= -^(  + e^  or  )H^  exp  i 

Using  the  Maxwell  equation  (6),  this  became 

(k^  + a^)  H. 

i € - 

The  parameter  was  thus  obtained! 


a = 


€ - k^ 


(56) 


(57) 


(58) 


Since  the  electric  field  was  expressed  in  terms  of  the 
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(59) 


magnetic  field,  the  electric  field  could  also  be  defined  as 


E = E^exp  i<l>+  + E_exp  i 4»_ 


(60) 


Thus  the  impedance 


A = 


®x  i «!>+  + E^_exp  i <p_ 


H. 


H^exp  i 4>_|_  1-  H_exp  i 9 


(61) 


At  the  upper  boundary,  '!>+  = 4>  = 0,  so  that 


_ aH^  - a H_ 
2=0  w€  (H^  + H^) 


(62) 


where  Eq  (57)  was  substituted  for  the  electric  field.  Thus, 
and  H_  are  related  by 


= 


+ A. 


-= A 

w€  *0  J 


H 


(63) 


where  A^  was  the  field  impedance  at  the  upper  boundary. 

Using  Eq  (57)  for  the  electric  field,  Eq  (63)  for  the  mag- 
netic field,  Eq  (59)  for  a , and  Eq  (62)  for  A^,  the  general 
expression  for  the  impedance  within  the  waveguide  was 


« (Aq  + + (Aq  - -^)exp  -i  a z 

A(2)  - (a  a )exp  iaz  + ( g - A )exp  -i  or  z 

L o o 


(64) 


Because  the  boundaries  were  infinite  conductors,  the  zero 
initial  condition  was  used  in  Eq  ( 64)  to  obtain  the  impedance 
at  the  lower  boundary. 


A(a)  = i(-jyp-)  tan  aa  = 0 


(65) 


where  a was  the  waveguide  depth.  Using  the  definition,  Eq  ( 33) t 
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this  suggested 


~T"  = wT-) 

(66) 

® z = u(  z) 

(67) 

For  the  impedance  to  be  zero, 

a a “ m-rr  = u(  a)  (68) 

Thus,  while  expression  of  the  boundary  condition  in  terms  of 
A produced  an  ambiguity  as  a result  of  the  complex  arg\iment 
of  the  tangent  function,  the  expression  in  terms  of  u gave 
the  boundary  condition  at  the  bottom  of  the  waveguide 
explicitly,  and  it  is  valid  even  for  bther  values  of' the 
constant  t , For  the  TM^^^  mode,  the  boundary  conditions 
for  Eq  (36)  become 

u(0)  = 0 and  u (a)  = TT  (69) 

The  expression,  Eq  (64),  also  could  be  used  to  check 
the  integration  by  the  computer.  Since  the  direct  integration 
solution  was  mathematically  identical  to  the  method  of 
deriving  Eq  (64),  the  values  obtained  by  each  for  a given 
k should  be  identical.  The  check  on  the  minimization  pro- 
cedure in  the  program  was  performed  using  the  value  of  k 
from  waveguide  theory  (Ref  22«185)» 

k = /Ac  - (f-)^  (70) 

This  was  compared  with  the  final  value  of  k produced  by  the 
computer  program. 
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£a?£-II.» Infinitely  Conducting  Boundaries  With  Small 

But  Finitely  Conducting  Waveguide  Medium.  The  primary 
difference  between  the  first  case  and  this  one  was  that  the 
dielectric  constant  of  the  waveguide  medium  became  complex: 


€ = 


(71) 


Thus  k became  complex^  However,  Eq  (4?),  the  upper  boundary 
condition,  remained  zero.  The  lower  boundary  condition 
remained  tt  , since  Eq  (68)  still  held.  Once  again  Eqs  (64) 
and  (70)  were  used  to  check  the  values  produced  by  the 
computer  program.  The  intemal  conductivity  was  defined  as 

<r  = e ' (72) 

Case  III:  Finitely  Conducting  Boundaries  with  Perfectly 
Dielectric  Waveguide  Medium.  The  initial  condition  for  this 
case  was  obtained  from  the  same  equation  as  that  for  the 
infinitely  conducting  boundaries: 

*o  ' (i  - 1)  (O?) 

However,  the  difference  was  that  since  the  conductivity  was 
not  infinite. 


6^0 


(73) 


Thus 


A 


o 


fa>  u 

2 <r 


(i 


1) 
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(74) 


Additionally,  A^  had  to  be  converted  to  the  phase  function 
form,  u^,  since  the  differential  equation  had  been  converted 
to  a function  of  u.  In  order  to  effect  the  conversion,  it 
became  necessary  to  solve  for  the  proportionality  constant. 
Since  i was  purely  imaginary,  a real  constant, 


X = -il 


was  defined.  Then  Eq  (36)  became 


du  - <«>€ 


{sin^  u + [1  - 


— ] cos^  u} 


(75) 


(76) 


X was  chosen  for  € real  and  k = k , such  that 

0 


^ C 1 - = 1 


(77) 


where 


ko  = - (J)" 


rn€’  - k. 


Thus  for  this  situation 


= const. 


z + u. 


(79) 


(81) 


which  was  the  linear  depth  dependence  desired.  The  phase 

kO 


function,  at  the  top  of  the  waveguide  (+),  or  at  the  bottom 
( -) , could  be  expressed 


u = tan”^( iAX) 


(82) 


= tan'^  C±‘Xyf^(-i  - 1)}  (83) 

= tan”^  [+X  exp  iTr/4}  (84) 

Using  a series  expansion  on  Eq  (82),  this  was  expressed 

u « iA'X  - (85) 


Factoring  out  ( iAx)  , and  using  Eq  (84), 

U = ±x  exp  ItrA  {1  . 11^  i . bii£pl  * tii!sx!jl2  1]  (86) 

50-'^  7<r'> 

The  function  P was  defined 


U W 
<r 


Also  the  expression 


exp  iw/4  = 


was  used.  Thus 


u = +p  [ 1 

V2 


1 i 

VF 


+ pli  ] 

3 5 7 


Therefore , if 


(87) 


(88) 


(89) 


u = X + iy 


(90) 


then  the  upper  boundary  conditions  were,  taking  the  (+)i-sign 
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in  Eqs  (83)  - (89) , 


P 

X = yf  {1  + y - 

p p^ 

y = ' 3 ■ 


(91) 

(92) 


Thus  expressions  (91).  (92),  and  (8?)  were  used  as  initial 
conditions. 

Assuming  a mode,  the  boundary  conditions  at  the 

lower  boundary  were,  taking  the  (-)  sign  in  Eqs  (83) -(89) 


P P^  P^  P^ 

* = ■ T " T “ 

p p^  p^  p^ 

y = -V5-C1  - y - y * - (94) 


The  -TT  difference  in  the  real  part  accounted  for  the  phase 
change  across  the  waveguide  for  the  lowest  mode.  It  should 
be  noted  that  P must  be  much  less  than  xmity  for  these  boundary 
conditions  to  hold.  For  certain  material  parameters,  it  would 
be  necessary  to  decrease  the  chosen  value  of  for  this 
expansion  to  be  valid. 

The  theoretical  k-values  for  a relatively  low-loss 
case  were  (Ref  22 t 184,  192) » 


/ <-3,2^ 

1 2r  ’ 

J 2<ra  Mif 

- (rr/a)^} 

- - (rr/a)^  (95) 


These  were  compared  with  those  obtained  from  the  computer 
calculation. 

Actual  Problem.  The  actual  problem  involved  the 
boundary  conditions  (91) -(9^)  with  a complex  conductivity. 
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The  conductivity  profiles  discussed  in  the  geophysics 
section  were  converted  to  functional  form  by  a curve-fit 
routine  on  the  computer,  A detailed  description  of  the 
process  has  been  placed  in  the  Appendix.  The  functional 
polynomial  then  served  as  the  imaginary  part  of  the  permit- 
tivity 

Im  (e)  = exp  (Az^  + Bz^  + Cz^  + Dz  + E)  (96) 

Because  no  previous  work  of  this  nature  has  been  found  on 
this  particular  problem,  no  real  comparison  of  the  final 
results  could  be  made  for  verification. 
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IV.  Noise  Analysis 


Introduction 

Surface  communications  must  deal  with  three  main 
noise  sources,  other  than  noise  inherent  to  transmitting  and 
receiving  devices:  atmospheric  noise,  produced  hy  electrical 
thunderstorms!  cosmic  noise,  produced  by  radio  stars  and 
solar  electrical  activity!  and  man-made  noise  produced  by 
other  transmitters  (Ref  3027).  The  primary  appeal  of 
thru-the -earth  communications  is  that  the  last  two  noise 
sources  can  be  totally  neglected  and  that  the  first  one  may 
be  attenuated  enough  to  be  virtually  neglected  (Ref  2:189). 

Thus,  there  are  only  two  external  noise  sources  in  the 
subsurface  waveguide:  atmospheric  noise  propagating  down 
through  the  crust  from  the  surface,  and  thermal  noise  due  to 
high  temperatures  within  the  crust.  The  original  motivation 
for  investigating  the  feasibility  of  lithospheric  communication 
was  the  hypothesis  that  increased  attenuation  of  transmitted 
signals  in  the  crust  might  be  more  than  compensated  for  by 
attenuation  of  the  atmospheric  noise  by  the  upper  litho- 
spheric strata,  to  quite  low  levels  (Ref  ktv) . A lower 
noise  level  than  at  the  surface  might  then  enable  the  sig- 
nals to  propagate  over  a significant  range.  This  hypothesis 
was,  of  course,  the  subject  of  this  study. 

The  basic  noise  problem  was  threefold:  to  calculate 
the  atmospheric  noise  field  in  the  suboceanic  waveguide  and 
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to  analyze  its  effects  on  the  reception  of  signals,  to  calcu- 
late the  thermal  noise  in  the  waveguide  and  its  effects  on 
reception,  and  to  calculate  the  signal-to -noise  ratio  and  the 
data  rate  possible  for  the  calculated  noise  levels.  These 
calculations,  along  with  the  propagation  constants  calculated 
from  waveguide  theory  in  Chapter  III , gave  an  estimate  of  the 
feasibility  of  using  the  suboceanic  crust  as  a communications 
channel . 

Atmospheric  Noise 

One  of  the  major  differences  between  the  subcontinental 
waveguide  and  the  suboceanic  waveguide  was  revealed  in  the 
analysis  of  the  attenuation  of  the  atmospheric  noise  as  it 
propagates  down  to  the  waveguide.  Although  the  overburden 
of  sedimentary  strata  on  the  continental  waveguide  serves  as 
a rather  effective  insulator  from  atmospheric  noise,  the 
combined  effects  of  several  kilometers  of  seawater  and  a 
shallow  sedimentary  overburden  on  the  oceanic  waveguide  are 
as  much  as  eight  times  the  insulator  as  the  continental 
overburden  (Ref  3*328). 

The  noise  density  spectrum  at  the  surface  has  been 
measured  by  E.L.  Maxwell  for  various  geographic  locations 
(see  Figure  7)  (Ref  23t64l).  Using  these  values  and  assumed 
constant  conductivity  values  for  various  materials  (seawater, 
granite,  sedimentary  rock),  Maxwell  also  calculated  the 
variation  of  the  noise  spectrum  with  depth  to  30O  meters. 
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Extrapolations  of  these  curves  (see  Figure  8)  produced  similar 

results  for  attenuation  values  at  1 tan  depth  as  those  given 

in  an  earlier  study  (see  Figure  9)  (Ref  3028).  Using  these 

extrapolations  to  a depth  of  4.7  km,  the  average  depth  of 

4 

the  oceans,  the  atmospheric  noise  at  a frequency  of  10  Hz, 
on  the  seafloor,  was  calculated  to  be  on  the  order  of 
times  the  magnitude  of  the  noise  at  the  surface.  Thus  it 
was  considered  negligible  at  the  ocean  floor,  and  even  more 
negligible  within  the  waveguide,  10-20  tan  below  the  surface 
of  the  crust.  Based  on  these  calculations,  the  atmospheric 
noise  was  totally  neglected  in  the  signal  attenuation  calcu- 
lations for  the  waveguide. 

Thermal  Noise 

The  second  source  of  noise  in  the  waveguide  is  due  to 
the  effects  of  high  temperatures.  The  relationship  between 
available  noise  power  and  the  temperature  of  a resistance 
in' a circuit  may  be  calculated  using  the  Norton  or  Thevenin 
equivalent  circuits  and  Nyquist's  Theoremi 

= kTB  (97) 

(Ref  24i469-470 ,481) . By  treating  the  waveguide  medium's 
thermal  properties  as  a resistor  in  a circuit,  the  thermal 
noise  may  be  calculated,  using  Eq  (97)  and  the  thermal 
profiles  discussed  in  the  Geophysical  Analysis  chapter. 

Detailed  analysis  and  calculations  of  the  noise  power 
have  been  performed,  which  take  into  account  the  directivity 
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P.D.  (dB  rel  to  1 A/m) 


Figure  9,  Depth  variation  of  horizontal  magnetic 

?i®ld  density  of  atmospheric  noise 

spectra  with  frequency.  Circles  are 
taken  from  Figure  2 (Ref  23i644),  and 
triangles  from  Ref  3*328,  Shows  the 
corroboration  of  the  two  sources. 
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of  the  antenna  and  thermal  contributions  from  the  surround- 
ing medium  (Ref  25»^-6) . However  these  detailed  results 
vary  only  by  one  or  two  decibels  from  the  results  of  Eq  (97) 

(Ref  4j25).  Thus  Eq  (97)  was  used  to  calculate  the  noise 
power  in  the  waveguide. 

The  thermal  noise  power  was  calculated  using  a thermal 
profile  deduced  from  heat-flow  measurements  performed  at  the 
seafloor  in  a high  heat-flow  region  (Ref  13*221).  The 
receiving  antenna  was  assumed  to  reach  to  a depth  of  15  km, 
where  it  encountered  a temperature  of  540°C.  Using  the  same 
bandwidth  as  in  the  atmospheric  noise  measurements,  6 Hz, 
the  noise  power  in  this  worst  case  was  calculated  to  be 
il62  dBm.  Since  the  thermal  noise  dominated  the  atmospheric 
noise,  the  thermal  noise  was  used  to  determine  signal-to- 
noise  ratios  at  the  receiving  antenna.  However,  as  often 
occurs  in  communications  systems,  the  thermal  noise  in  the 
receiver  was  considered  to  dominate  even  the  thermal  noise 
in  the  waveguide,  and  the  sum  was  arbitrarily  set  at  -160  dBm. 

This  is  still  a significant  improvement  over  the  noise  situ- 
ation in  surface  communications,  where  receiver  thermal 
noise  is  often  dominated  by  atmospherics  and  external  thermal 
noise  (Ref  3*328). 

System  Information  Capacity 

The  details  of  pulse-code  modulation  (PCM)  of  a signal 
in  a communications  system  will  be  discussed  here.  Interested 
readers  are  referred  to  the  excellent  description  in 
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Ref  26*138-163.  It  is  sufficient  to  state  that  the  information 
capacity  of  a system  is  dependent  on  the  minimum  response 
time  of  the  system  and  the  number  of  distinguishable  voltage 
levels  in  the  coding  scheme.  Using  the  average  power  over 
a large  interval,  the  maxiumum  rate  of  information  trans- 
mission is 

C = B logg  (1  + |)  (98) 

where  S/^  is  the  signal-to-noise  ratio  (SNR)  for  the  system 
(Ref  26*150).  The  required  SNR  cannot  be  determined  without 
assuming  certain  factors  about  the  communication  devices. 

Given  a desired  bandwidth  and  data  rate  capacity  for  the 
above  equation,  the  required  SNR  can  be  calculated.  For  a 
rate  of  one  bit/second  and  bandwidth  of  one  hertz,  the  SNR 
should  be  at  least  unity.  In  practice,  the  SNR  often  must 
be  several  dB  above  the  limit  given  in  Eq  (98)  (Ref  25*24). 

The  results  of  the  waveguide  calculations,  the  noise 
power,  and  the  information  capacity  can  be  used  to  calculate 
the  necessary  transmitted  power  to  obtain  an  assumed  data 
transmission  capacity.  Since  the  total  noise  power  is  the 
sum  of  the  thermal  and  receiver  noise,  let 

N = N^  + N^  (99) 

and  the  receiver  signal  power  can  be  related  to  the  trans- 
mitted power  by 


S = S^  Cexp  '2  Im(lac)} 
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r 

whera  is  the  transmitted  power.  The  required  signal 
power  at  the  transmitter  is 

Sq  = Cexp  2 Im(kx)}  (N^  + - 1)  (101) 

for  a given  information  rate  and  bandwidth.  In  this  wayi  it 
was  at  least  theoretically. possible  to  calculate  the  feasi- 
bility of  transmitting  a signal  of  sufficient  power  to  be 
received  and  decoded  after  travelling  a desired  range,  x. 
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V . Data  and  Results 


Input  Data 

Performance  of  the  calculations  for  the  attenuation 
coefficient  required  the  use  of  two  fundamental  values:  the 
magnetic  permeability  , and  the  real  electric  permittivity  of 
the  waveguide  medium.  The  magnetic  permeability  was  assumed 
to  be  identical  to  that  of  free  space,  since  the  materials 
of  not  only  the  waveguide  medium  but  also  of  the  boundaries 
were  assumed  non-magnetic  (Ref  32:2).  The  value  of  the 
electric  permittivity  used  in  calculation  varied  in  the 
literature.  The  value  used  in  this  study  was  farad/ 

meter  (Ref  10:5-137) » since  the  suboceanic  waveguide  medium 
seemed  to  be  composed  primarily  of  basalt.  The  depth  of  the 
waveguide  varied  from  one  conductivity  profile  to  another, 
due  to  the  nature  of  the  computer  routine  LOOKAT.  This  pro- 
gram stopped  calculations  at  the  depth  when  the  conductivity 
had  increased  from  its  minimum  value  to  one  where  the 
imaginary  part  of  the  dielectric  constant  was  within  an 
order  of  magnitude  of  the  real  part.  In  other  words,  the 
calculation  stopped  when  conductor  losses  became  large  enough 
to  prevent  propagation  in  the  waveguide.  For  the  shallow 
oceanic  profile,  this  occurred  irt  the  vicinity  of  30  km. 

The  conductivity  profiles  were  input  into  the  program 
in  the  form  of  a polynomial,  Eq  (96).  The  coefficients  of 
this  polynomial,  calculated  by  the  program  FIT,  are  listed 
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in  Table  I,  with  further  details  given  in  the  Appendix. 

Output  Data 

Because  of  severe  time  constraints,  only  the  shallow 
oceanic  profile  was  treated  fully  for  the  TM^q  mode.  Calcu- 
lation of  the  final  propagation  and  attenuation  constants 
for  one  wavelength  value  involved  perfoming  successive 
computer  runs,  adjusting  the  upper  and  lower  limits  on  k, 
until  the  values  began  to  converge,  with  as  small  an  error 
from  the  boundary  conditions  as  possible.  Then  a final  run 
was  made,  using  a large  number  of  convergent  steps,  to 
calculate  a final  value  of  k.  The  values  calculated  for  the 
shallow  oceanic  profile  are  listed  in  Table  II.  These  values 
were  then  used  to  calculate  the  attenuation  within  the  wave- 
guide, when  used  as  a communications  system. 

Results 

Using  a quadratic  curve-fit  procedure,  described  in 
the  Appendix,  the  minimum  attenuation  was  calculated  to  be 
1»55^7  X 10“^  nepers/meter  at  a frequency  of  1.3  x 10^  Hz. 
Assuming  a transmitted  power  of  1 kW  for  the  planar  wave 

Table  II.  The  propagation  constant,  k = 3 + ia,  for  various 
frequencies,  and  the  respective  error  in  the 
calculation. 


(Hz) 

B (n/m) 

g (n/m) 

Error(  rel) 

.9.0 

j 

6.3787  X 10"^ 

5.7673  X 10’^ 

5.6215  X 10'^ 

1.1  X 10^ 

8.4795  X 10~5 

2.4229  X 10'^ 

4.8118  X 10"^ 

1.5  X 10^ 

6.4953  X 10"^ 

3.4329  X 10“^ 

7.4953  X 10"^ 

1.9  X 10^ 

1.3162  X 10"^ 

4.5394  X 10“^ 

7.7851  X 10"^ 
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Coefficients  for  conductivity  polynomial,  calculated 
in  FIT,  for  the  different  conductivity  profiles. 


Table  I. 

Mott  & Biggs  ( deep) 

A = -7.8826  X 10"^® 
B = 4.4923  X 10"^^ 

C = -3.1180  X 10"^ 

D = 6.3691  X 10"^ 

E = -1.6087  X 10^ 

Hauselv  (EUS) 

A = 2.6085  X lO"^"^ 

B = -3.5396  X 10'^^ 
C = 1.7246  X lO""^ 

D = -3.2968  X 10'^ 

E = 1.5213  X 10“^ 

Levin  ( dry) 

A = 3.1109  X 10"^^ 

B = -2.8981  X 10“^^ 
C - 1.1069  X lO'"^ 

D = -1.7696  X 10"^ 

E = -5.8058  X 10° 


Mott  & Biggs  ( shallow) 

A = -9.6638  X 10"^® 

B = 7.O8O5  X 10"^^ 

C = -1.2495  X 10"® 

D = 8.4344  X 10"^ 

E = -1.3861  X 10^ 

Houselv  (BR) 

A = 8.1033  X 10’^^ 

B = -7.4797  X 10'^^ 

C = 2.7101  X lO""^ 

D = -4.0537  X 10'^ 

E = I.8II9  X 10° 

Levin  (wet) 

A = 9.4523  X lO"^"^ 

B = -6.6532  X 10'^^ 

c = 1.5513  X lO"”^ 

D = -1.4071  X 10"^ 

E = -6.7000  X 10° 


Keller 

A = 5.6825  X lO"-*-^ 

B = 9.8938  X 10"^^ 

G = -1.2703  X lO""^ 

D = -5.0319  X 10’^ 
E = -9.2258  X 10° 
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order  to  obtain  a minimum  SNR  of  unity.  The  minimum 
attenuation,  1.55^7  x 10”^  nepers/meter,  was  calculated  to 
be  -1.35  dB  /Van.  Therefore  the  meiximum  range  of  the  system, 
under  these  circumstances,  is  I63  Icilometers. 

Attempts  to  calculate  the  attenuation  using  a continental 
profile  were  unsuccessful.  Because  of  time  constrainst  on 
the  study,  the  errors  preventing  this  set  of  calculations,  as 
well  as  those  for  the  quasi-TM^Q  mode  were  not  eliminated. 


VI . Conclusions  and  Recommendations 


TKIs  report  has  attempted  to  summarize  the  work 

perfomed  in  a study  of  the  earth's  crust  under  the  ocean  as 

a very  low  frequency  communications  channel.  The  study  was 

broken  into  three  parts;  a geophysical,  waveguide^  and  noise 

analysis.  The  goephysical  analysis  examined  the  effects 

of  structural,  thermal,  and  electrical  properties  on  the 

conductivity-depth  profile  of  the  crust.  This  profile  was 

then  applied  to  an  analysis  of  the  electromagnetics  of  the 

crust  as  a waveguide , to  derive  and  calculate  the  attenuation 

\ 

constant  for  a planar  mode.  The  atmospheric  noise  was 

neglected  due  to  its  attenuation  by  the  overlying  ocean,  and 
the  thermal  noise  due  to  temperature  effects  in  the  waveguide 
was  calculated  to  be  -162  dBm.  The  combined  receiver  and 
thermal  noise  was  chosen  to  be  -l6o  dBm. 

The  attenuation  in  the  waveguide  was  calculated  to  be 
-1.35  dB  /km,  which  produced  a maximum  range  of  about  I63  km 
for  a minimum  signal-to-no3 se  ratio  of  unity.  It  should  be 
noted  that  these  results  neglect  receiver  noise,  and  there- 
fore represent  an  upper  limit  on  range.  The  attenuation 
found  here  was  corroborated  by  several  other  sources  (Refs 
3*326,  4:33,  32*54),  to  within  an  order  of  magnitude.  On 
the  basis  of  these  results,  it  is  not  feasible  to  think  in 
terms  of  using  the  earth’s  crust  for  a very  low  frequency 
communications  system. 
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Although  the  limitations  to  this  study  could  be  used 
to  justify  further  feasibility  studies  on  this  subject,  this 
is  not  recommended.  In  fact,  there  have  already  been  enough 
studies  of  the  various  aspects  of  the  lithospheric  commimi- 
cations  problem  from  the  waveguide  perspective  to  justify 
discarding  the  idea  altogether.  The  reason  the  idea  has 
persisted  over  the  years  since  its  introduction  is  its 
aesthetic  appeal,  as  well  as  the  fact  that  it  fills  a very 
real  need.  However,  this  system  may  well  be  compared  to  the 
early  trans-Atlantic  cables — because  of  technical  problems, 
the  initial  concept  of  simply  a telegraph  wire  across  the 
Atlantic  did  not  work.  The  system  simply  was  not  feasible 
until  the  later  innovations  in  insulation  and  repeater 
technology  came  about.  Therefore,  rather  than  persisting 
in  trying  to  make  a lithospheric  waveguide  communications 
system  work  by  "brute  force"  as  in  this  study,  it  would  be 
better  to  start  with  the  basic  idea  of  a hardened  system, 
and  develop  it  along  more  innovative  lines.  There  may 
indeed  be  a way  to  use  the  earth’s  crust  for  communications, 
but  the  indications  are  that  a very  low  frequency  waveguide 
is  not  one  of  them. 


58 


1 


1. 


2. 


3. 


5. 


6. 


7. 


8. 


9. 

10. 


Bibliography 


Hibben,  Stuart  G.  "Underground  Electromametic  Wave 
Propagation,"  Foreign  Science  Bulletin,  ^(3) *20-32 
(March  196?) • 

Wheeler,  Harold  A.,  "Radio  Wave  Propagation  in  the 
Earth ’ s Crust , " Journal  of  Research  of  the  National 
Bureau  of  Standards — D.  Radio  Propagation,~S5D(  2)7  189- 
191  ( March-April  19617, 

Mott,  H.  and  A.W.  Biggs,  "Very-Low -Frequency  Propagation 
Below  the  Bottom  of  the  Sea,"  IEEE  Transactions  on 
Antennas  and  Propagation . AP-11; 323-329  (May  19637. 

Field,  E.C.  and  M.Dore,  Electromagnetic  Communication 
in  the  Earth ' s Crust . Pacific-Sierra  Research  Corp. , 
Report  306,  Advanced  Research  Projects  Agency,  ARPA 
Order  No.  2270,  June  1973.  (AD  913  801) 

Chadwick,  P.,  "Mountain  Building  Hypotheses,"  Interna- 
tional Geophysics  Series . Volume  2,  Continental  Drift , 
edited  by  S.K.  Runcorn.  New  York  and  London*  Academic 
Press,  1962. 

Spencer,  Edgar  W.  Introduction  to  the  Structure  of  the 
Earth . New  York:  McGraw-Hil.L  Books,  I969 . 

Heezen,  Bruce  C.,  "The  Deep-Sea  Floor,"  International 
Geophysics  Series , Volume  Continental  Drift , edited 
by  S.K.  Runcorn.  New  York  and  London:  Academic  Press, 

1962. 

Zvereva,  Ye.V.,  A.M.  Ryazantsev,  I .N . Samuylov,  and  D. 

N . Shakhsuvarov . Research  on  the  Propagation  of  Electro- 
magnetic Waves  in  Earth  * s Crust.  Working  Copy: 

Machine  Translation.  Wright-Patterson  AFB,  Ohio: 

Foreign  Technology  Division,  25  February  1976.  (FTD- 
ID(RS)T-0183-76) . Originally  published  as  "Izvcheniye 
Rasprostraneniye  Yelektromagnitnykh  Voln  v Zemnoy  Kore,"  , 
Rasprostrainenive  Radiovoln . Moskva:  Izdatel'stvo 
"Nauka" , 1975.  Both  original  and  translation  used  by 
the  author. 

Weast,  Robert  C.,  editor.  Handbook  of  Chemistry  and 
Physics  (fifty-second  edition) . ^eveland,  Ohio: 

Chemical  Rubber  Company,  1971. 

Gray,  Dwight  E.,  coordinating  editor.  American  Institute 
of  Physics  Handbook  (third  edition).  New  York*  McGraw- 


59 


Hill  Book  Company,  1972. 


11.  Nelson,  A.  and  K.D.  Nelson.  Dictionary  of  Applied 
Geology,  Mining,  and  Ciyil  Engineering.  New  York: 
Philosophical  Library,  Inc.,  I967. 

12.  Cox,  Charles  S.,"The  Electrical  Conduct iyity  of  the 
Ocesuiic  Lithosphere,"  Geophysical  Monograph  Series . 

Volume  14,  The  Structure  and  Physical  Properties  of  the 
Earth ' s Crust,  edited  by  John  G.  Heacock.  Washington, 
D.C.i  American  Geophysical  Union,  1971. 

13.  von  Herzen,  R.P.,  "Surface  Heat  Flow  and  Some  Implica- 
tions for  the  Mantle,"  The  Earth  * s Mantle . edited  by 
Thomas  F.  Gaskell.  London  and  New  York:  Academic 
Press,  1967. 

14.  Lubimoya,  E.A. , "Theory  of  Thermal  State  of  the  Earth's 
Mantle,"  The  Earth's  Mantle . edited  by  Thomas  F.  Gaskell. 
London  and  New  York:  Academic  Press,  1967- 

15.  Beck,  Alan  E ., "Techniques  of  Measuring  Heat  Flow  on 
Land,"  Geophysical  Monograph  Series . Number  8, 
Terrestrial  Heat  Flow,  edited  by  William  H.K.  Lee. 
Washington,  D.C.:  American  Geophysical  Union  of  the 
National  Academy  of  Sciences--National  Research  Council, 

1965. 

16.  Blackwell,  Dayid  D.,  "The  Thermal  Structure  of  the 
Continental  Crust,"  Geophysical  Monograph  Series , Volume 
l4.  The  Structure  and  Physical  Properties  of  the  Earth's 
Crust,  edited  by  John  G.  Heacock.  Washington,  D.C.: 
American  Geophysical  Union  1971. 

17.  Makaroy,  G.I.,  and  V.  A.  Payloy,  "Suryey  of  Works 
Connected  with  Underground  Radiowaye  Propagation," 
Unedited  Machine  Translation.  Wright-Patterson  AFB, 

Ohio:  Foreign  Technology  Diyision,  27  July  1977. 

(FTD-ID( RS)T-1243-77)  Originally  published  as  "Obzor 
Rabot,  Syyazanykh  s Podzemnym  Rasprostraneniyem  Radio- 
voln,"  Rasprostraneniye  Radioyoln,  Vypusk  IV,  Problemy 
Difraktsii  i,  Rasprostraneniya  Voln.  Leningrad:  Lenin- 
gradskogo  Uniyersiteta,  i960.  Both  original  and  trans- 
lation used  by  the  author. 

18.  Keller,  G.V.,  "Electrical  Properties  in  the  Deep  Crust," 
IEEE  Transactions  on  Antennas  and  Propagation , AP-11 : 
344-357  (May  1963) . 

19.  Levin,  S.  Benedict,  "Lithospheric  Radio  Propagation — 

A Review,"  A CARD  Conference  Proceedings  No . 20 , Selected 
Papers  from  the  12th  Symposium  of  the  Ayionics  Panel 


60 


I 


of  AGARD,  held  in  Paris,  France,  25-29  April  I966. 

Londoni  Technical  Editing  and  Reproduction,  Ltd., 

1966. 

20.  Garland,  G.D.  and  T.F.  Webster,  "Studies  of  Natural 
Electric  and  Magnetic  Fields,"  Journal  of  Research  of 
the  National  Bureau  of  Standards — D.  Radio  Propagation. 
^4D(  4m  40  5-^08  (Julv-August  IQ60)  . 

21.  Levin,  S .Benedict , "Radio  Propagation  Through  the  Crust — 
Retrospect  and  Prospect,"  Geophysical  Monograph  Series . 
Volume  14,  The  Structure  and  Physical  Properties  of  the 
Earth's  Crust . edited  by  John  G.  Heacock.  Wasington, 

D.C.i  American  Geophysical  Union,  1971. 

22.  Jordan,  Edward  C.  and  Keith  G.  Balmain.  Electromagnetic 
Waves  and  Radiating  Systems  ( second  editionT"!  Engle- 
wood  Cliffs,  New  Jersey:  Prentice-Hall,  Inc.,  I908. 

23.  Maxwell,  Eugene  L.,  "Atmospheric  Noise  From  20  Hz  to 

30  Idiz,"  Radio  Science  (New  Series),  2(6):637-644  (June 

1967). 

24.  Ziemer,  R.E.  and  W.H.  Tranter,  Principles  of  Communi- 
cations; Systems . Modulation,  and  Noise . Boston: 

Houghton  Mifflin  Co . , I976 . 

25.  Gallawa,  R.L.  and  L.L.  Haidle.  ^ Engineering  Feasi- 
bility Study  of  Communications  Via  a Subterranean 
Waveguide . Telecommunications  Research  and  Engineering 
Report  39*  Boulder  Colorado:  U.S.  Department  of 
Commerce,  Office  of  Telecommunications,  Institute  for 
Telecommunications  Sciences,  November  1972.  (AD  75^  786) 

26.  Schwartz,  Mischa.  Information . Transmission.  Modulation. 
pid  Noise . A Unified  Approach  to  Communications  Systems 

( second  edition) .New  York:  McGraw-Hill  Book  Co.,  1970. 

27.  Garland,  G.D.  and  T.F.  Webster,  "Studies  of  Natural 
Electric  and  Magnetic  Fields,"  Proceedings  of  i960 
Conference  on  the  Propagation  of  ELF  Radio  Waves. 

Boulder,  Colorado:  National  Bureau  of  Standards, 

January  i960. 

28.  Hughes,  H.,  "The  Electrical  Conductivity  of  the  Earth's 
Interior,"  Ph.D.  dissertation.  University  of  Cambridge, 
Cambridge,  England,  1953. 

29.  Simmons,  G. , R.M.  Housely,  and  F.J.  Forin,  "Electrical 
Resistivity  of  Rocks  and  the  Possibility  of  Lithospheric 
Communication,"  to  have  been  published  in  1973.  (Private 


61 


communication  with  authors  indicates  that  it  was  never 
published. ) 

30.  Schwering,  F.K.,  D.W.  Peterson,  and  S.B.  Levin,  "A 

Model  for  Electromagnetic  Propagation  in  the  Lithosphere," 
Proceedings  of  the  IEEE , _56:799-804  (May  I968). 

31.  Brace,  W.F.,  "Resistivity  of  Saturated  Crustal  Rocks 

to  40  km  Based  on  Laboratory  Measurements,"  Geophysical 
Monograph  Series , Volume  14.  The  Structure  and  Physical 
Properties  of  the  Earth ' s Crust , edited  by  John  G. 

Heacock.  Washington,  D.C.:  American  Geophysical 
Union,  1971. 

32.  Spies,  Kenneth  P.  and  James  R.  Wait.  On  Calculations 
of  the  Modal  Parameters  of  ^ Idealized  Earth-Crust 
Waveguide.  Technical  Report  No.  1,  Contract  No.  NA-ONR- 
15-71.  Boulder,  Colorado:  U.S.  Department  of  Commerce, 
Office  of  Telecommunications,  Institute  for  Telecommuni- 
cations Sciences,  17  March  1971. 


Supplementary  Bibliography 


General 

Banning,  C . J . , K.R.  Cook,  and  R.D.  Kelly,  Experimental  and 
Theoretical  Inyestigations  of  the  Feasibility  of 
Waye  Propagation  in  the  Earth ' s Crust . Technical 
Report  EE-83.  Albequerque,  New  Mexico s The  Uniyer- 
sity  of  New  Mexico,  October  1962.  (AD  290  776) 

Bays,  C.A.  Study  of  the  Propagation  of  Electromagnetic  Wayes 
Through  Lithologic  Formations . Final  Progress  Report, 
Contract  No.  DA  36 -039 -SC -5^5^.  Urbana,  Illinois: 

C.A.  Bays  and  Associates,  Inc.,  195^.  (AD  33  576) 

Briscoe,  H.W.  Study  of  the  Feasibility  of  Lithospheric 
Electromagnetic  Communications , Final  Report, 

Report  No.  2346.  Bolt,  Beranek,  and  Newman,  Inc., 
Contract  No.  N00014-72-C-0310 , I6  May  1972.  (AD  894 
807) 

Burgess,  B.,  "Propagation  of  VLF  Wayes  Under  Disturbed  Condi- 
tions," Radio  Science  J oumal  of  Research,  NBS/\JSNC- 
URSI,  68D( 1) : January  1964. 

Burrows,  C.R.,  "Radio  Communication  Within  the  Earth's  Crust," 
IEEE  Transactions  on  Antennas  and  Propagation . AP-11 
(377311-317  (May  I963TT 

Cooper,  H.F.,  M.K.  Groyer,  and  D.B.  Wilson.  Description  and 

Possible  Capability  of  a Superhard  Underground  Command 
System  (SHUCS)  and  Associated  Technology  Needs , 

Topical  Report.  Marina  Del  Rey,  California:  R & D 
Associates,  Contract  No.  DNA  001-75-C-0005»  June 
1975.  (AD  BOll  654) 

Crain,  C.M.  An  Oyeryjew  Discussion  of  Propagation  Effects  of 

Nuclear  Enyironments  on  VLF -LF  Communications  Systems , 
Topical  Report  for  Period  4 Sept  74-31  Aug  75.  Santa 
Monica,  California:  The  Rand  Corporation,  Contract 
Na  DNA  OOI-74-C-OOI9,  31  August  1975.  (AD  A024  956) 

deBettencourt , J.T.  Review  of  Radio  Propagation  Below  the 
Earth ' s Surface . Paper  presented  at  XVth  General 
Assembly  of  URSI , Munich,  Germany,  August  I966. 

Prieraan,  E.A.  and  N.M.  Kroll.  Lithospheric  Propagation  for 
Undersea  Communication.  Technical  Report  JSR-73-5. 

Menlo  Park^  California:  Stanford  Research  Institute, 
Contract  No.  DAHC15-73-C-0370 , November  1973. (AD  917  802) 


63 


I 


I 


L 


Gallawa,  R.L.  and  L.L.  Haidle.  An  Engineering  Feasibility 

Study  o£_  Communications  yia  a Subterranean  Wayeguide. 
tRER-39 . Boulder,  Colorado:  Institute  for  Telecom- 
munications Sciences,  U.S.  Department  of  Commerce, 

November  1972.  (AD  75^  786) 

Hales,  A.L.  Crustal  Studies  Workshop:  A Multidisciplinary 
Program  of  Investigation  of  the  Possibility  of  an 
Electromagnetic  Waveguide  in  the  Lower  Crust . Dallas, 

Texas:  University  of  Texas,  Contract  No.  N00014-67- 
A-0310-0006,  7 March  1972.  (AD  7^3  1^6) 

Levin,  S. Benedict.  Lithospheric  Radio  Propagation- -A  Review. 

Fort  Monmouth,  New  Jersey:  Institute  for  Exploratory 
Research,  U.S.  Army  Electronics  Command,  i960  (?). 

(AD  701  017) 

O'Neal,  I.P.,  M.J.  McGowan,  F.  Kelly,  and  F.J.  Rhoads,  VLF 
Propagation  Predictions  for  Disturbed  Conditions , 

Report  8017,  Naval  Research  Labs,  Electromagnetic 
Propagation  Division,  Communication  Sciences  Branch, 

27  July  1976. 

Petrovskiy,  A.D.,  "Formulating  the  Theory  of  Underground 

Radio  Wave  Methods,"  Prikladnava  Geof izike , No.  57* 

155-185,  in  Russian. 

Ryazantsov,  A.M.  and  A.  Shabelnikov,  "Propagation  of  Radio 
Waves  Through  the  Earth's  Crust  (A  Review),"  Radio- 
technika  i Electronika  (Moscow),  10:1923  ( 195^TT 

Spies,  Kenneth  P.  and  J.R.  Wait.  Attenuation  Calculations 

for  a Two -Layer  Earth-Crust  Wave-Guide  for  Sub-Surface 
Electromagnetic  Propagation . Technical  Report  No. 2. 

Boulder,  Colorado:  Institute  for  Telecommunications 
Sciences,  Office  of  Telecommunications,  U.S.  Department 
of  Commerce,  Contract  No.  NA-ONR-15-71 • (AD  726  75^) 

Walt,  James  R.,  "An  Anomalous  Propagation  of  Radio  Waves  in 

Earth  Strata,"  Geophysics.  12.(2):  1954.  i 

, "Note  on  Calculations  of  Propagation  Parameters  for 

an  Idealized  Earth-Crust  Waveguide,"  Geophysical 
Mon graph  Series , Volume  14.  The  Structure  and  Physical 
Properties  of  the  Earth's  Crust , edited  by  John  G. 

Heacock.  Washington,  D.C.:  American  Geophysical 
Union,  1971. 

• On  Mode  Conversion  in  the  Earth-Crust  Waveguide , 

Technical  Report  No.~^.  Boulder,  Colorado:  Cooper- 
ative Institute  for  Research  in  Environmental  Sciences, 

University  of  Colorado,  and  Institute  for  Telecommun- 

64 


ication  Sciences,  Office  of  Telecommunications,  U.S. 
Department  of  Commerce,  Contract  No.  NA-ONR-15-71 , 

27  October  1972.  (AD  756  693) 


, "The  Possibility  of  Guided  Electromagnetic  Waves  in 
the  Earth's  Crust,"  IEEE  Transactions  on  Antennas 
and  Propagation.  AP-lK  3)  t 330-33*^  (May  1963)  • 

, "Propagation  of  Electromagnetic  Pulses  in  Terrestrial 
Waveguides,"  IEEE  Transations  on  Antennas  and  Propa- 
gation , AP-13(6);  901-918  (November  1965).  (AD  628 

Wo) 

, "Some  Factors  Concerning  Electromagnetic  Wave  Propa- 
gation in  the  Earth’s  Crust,"  Proceedings  of  the 
IEEE.  ^(8):  1020  ( 1966). 

,"A  Summary  of  VLF  and  ELF  Propagation  Research," 
Journal  of  Research  of  the  National  Bureau  of 
Standards . D.  Radio  Propagation.  5^D(  6) ; 647^^9 
(December  I960TT 

, "Terrestrial  Propagation  of  VLF  Radio  Waves  (A 
Theoretical  Investigation),"  Journal  of  Research  of 
the  National  Bureau  of  Standards . D.  Radio  Propagation . 
^(2)  j153-203(  March-April  i960). 

• Underground  Electromagnetic  Propagation,  Preliminary 
Results  for  the  period  1 Aug  70-31  Dec  70,  ARPA  Order 
No.  1361.  Boulder,  Colorado;  NOAA  Environmental 
Research  Laboratories;  and  Office  of  Telecommunica- 
tions, Institute  for  Telecommunications  Sciences, 
Contract  No.  910-1468,  15  January  1971.  (AD  882  757) 


Geophysics 

Ames,  L.A.,  J.W.  Frazier,  and  A.S.  Orange , "Geological  and 
Geophysical  Considerations  in  Radio  Propagation 
Through  the  Earth's  Crust,"  IEEE  Transactions  on 
Antennas  and  Propagation.  AP-11(  3) 1 369-371  (May  I963) . 

Beitchman,  J.G.  Soviet  Micro pulsations  Research  for  Possible 
Communications  Application.  Advanced  Projects 
Research  Administration,  ARPA  No.  2225,  21  February 
1973.  (AD  759  380). 

Berdichevskiy,  M.,  L.L.  Van’yan,  and  E.B.  Faynberg,  "Theo- 
retical Principles  in  Using  Electromagnetic  Variations 
to  Study  the  Electric -Conductivity  of  the  Earth," 
Geomagnetism  and  Aeronomy . IX(  3) 1 46 5-467. 

65 


L 


Bitterlich,  Wolfgang.  Measurement  of  the  Conductivity  and 

Dielectric  Constant  of  a Homogeneous  Medium  by  Means 
of  Wave  Propagation . Scientific  Report  No.  7i  AFCRL- 
. Innsbruck,  Austria:  VLF  Project,  2$  July 
1968.  (AD  677  587) 

Bloomquist,  M.,  H.W.  Smith,  F.X.  Bostick,  , Propagation 

Characteristics  of  Micropulsations , Final  Report, 
AFCRL-67-05^0 . Austin,  Texas:  Electrical  Engineer- 
ing Research  Laboratory,  University  of  Texas,  Contract 
No.  19(628)-4o46,  1 August  I967.  (AD  66I  284) 

Boettcher,  A.L,  "The  Nature  of  the  Crust  of  the  Earth,  with 
Special  Emphasis  on  the  role  of  Plagioclase , " Geo- 
physical Monograph  Series.  Volume  l4.  The  Structure 
and  Physical  Properties  of  the  Earth ' s Crust,  edited 
by  John  G.  Heaccc.v.  Washington,  D.C.:  American 
Geophysical  Union,  1971. 

Bondarchuk,  V.G.  Structure  of  the  Earth ' s Crust . Kiev: 

Publishing  House  of  the  Academy  of  Science  of  the 
Uktainian  SSR,  I962.  Translated  by  George  Synenko , 

TC  846.  St.  Louis,  Missouri:  Linguistic  Section, 
Aeronautical  Chart  and  Information  Center,  April 

1965.  (AD  623  443) 

Bostick,  F.X.,  D.R.  Word,  and  H.W.  Smith,  "Crustal  Investi- 
gations by  the  Magnetotelluric  Tensor  Impedance 
Method,"  Geophysical  Monograph  Series . Volume  14. 

The  Structure  and  Physical  Properties  of  the  Earth  * s 
Crust . edited  by  John  G.  Heacock.  Washington,  D.C.; 
American  Geophysical  Union,  1971. 

Braile,  L.W.  Investigation  of  Seismic  Velocity  Transition 
Zones  and  Application  to  the  Moho  in  the  Western 
United  States.  West  Lafayette,  Indiana:  Purdue 
University,  Department  of  Geosciences,  Contract  No. 
NOOOI4-75-C-O972,  30  June  1975.  (AD  012  542) 

Dedkova,  D.,  J.  Halousek,  B.  Krcmar,  and  K.  Prihoda,  "Geo- 
thermal Prospecting  in  Shallow  Holes  and  its  Limita- 
tions," Geothermics . Special  Issue  2,  U.N.  Symposium 
on  the  Development  and  Utilization  of  Geothermal 
Resources,  2:1244-1249  ( 1970) . 

Forsyth,  D.W.  and  F.  Press,  "Geophysical  Tests  of  Petrological 
Models  of  the  Spreading  Lithosphere,"  Journal  of 
Geophysical  Research.  76(  32) 1 7963-7979  (November 
10,  1971).  (AD  739  65^ 

Gangi,  A.F.  and  N.E.  Lamping, "An  in  Situ  Method  of  Determin- 


66 


ing  the  Pressure  Dependence  of  Phase-Transition 
Temperatures  in  the  Crust,"  Geophysical  Monograph 
Series . Volume  14.  The  Structure  and  Physical  Prop- 
erties of  the  Earth's  Crust,  edited  by  John  G. 

Heacock.  Washington,  D.C.:  American  Geophysical 
Union,  1971* 

Heacock,  John  G. , editor.  The  Structure  and  Physical  Prop- 
erties of  the  Earth ' s Crust , Geophysical  Monograph 
Series . Volume~l4.  Washington,  D.C.:  American 
Geophysical  Union, 1971. 

Kaspar,  M.  and  J.  Pecen,  "Detection  of  Cayes  in  a Karst  For- 
mation by  Means  of  Electromagnetic  Waves,"  Geophysical 
Prospecting.  23( 4) i 611-621 . 

Keller,  G.V. , "Electrical  Studies  of  the  Crust  and  Upper 
Mantle,"  Geophysical  Monograph  Series . Volume  l4, 

The  Structure  and  Physical  Properties  of  the  Earth ' s 
Crust,  edited  by  John  G.  Heacock.  Washington,  D.C.: 
American  Geophysical  Union,  1971- 

Keller,  G.V.  Workshop--Detennination  of  the  Minimum  Conduc- 
tivity of  the  Earth ' s Crust , Final  Report.  Golden, 
Colorado:  Colorado  School  of  Mines,  Contract  No. 
N00014-71-A-0430-0002,  23  June  1972.  (AD  743  867) 

Kellner,  W.  Electrical  Rock  Characteristics  From  the  Field 

Structure  of  a Magnetic  Dipole  Immersed  in  a Conduct- 
ing Medium.  Scientific  Report  No.  1,  AFCRL-S9-0281 . 
Innsbruck,  Austria;  VLF  Project,  Contract  No.  F 
6IO52-69-C-OOO7,  15  April  1969.  (AD  690  167) 

Khalafalla,  A.  Dielectric  Relaxation  and  Rock  Geophysical 
Characteristics . Semi-Annual  Report.  St.  Paul, 
Minnesota:  Research  Department,  Systems  and  Research 
Center,  Honeywell,  Inc.,  Contract  No.  H0210026, 

15  September  1971*  (AD  732  024) 

Knopoff , L.  and  Frank  Press.  Study  of  the  Earth ' s Crust  in 
Selected  Regions  of  the  World . Final  Report,  AFOSR 
4803 • Los  Angeles,  California:  Institute  of  Geo- 
physics, University  of  California;  Pasedena,  Cali- 
fornia: Seismological  Laboratory,  California  Institute 
of  Technology,  15  September  I962.  (AD  407  054) 

Madden,  T.R.,  "The  Resolving  Power  of  Geoelectric  Measure- 
ments for  Delineating  Resistive  Zones  within  the 
Crust,"  Geophysical  Monograph  Series . Volume  l4. 

The  Structure  and  Physical  Properties  of  the  Earth ' s 
Crust,  edited  by  John  G.  Heacock.  Washington,  D.C.: 
American  Geophysical  Union,  1971 • 


67 


Nessler,  N.,  M.  Oberladstatter , and  Dieter  Fritz.  Frequency 
Dependence  of  Electrical  Rock  Parameters  in  the  VLF 
Range . Final  Report.  Innsbruck,  Austria:  Electronics 
Laboratory,  Innsbruck  University,  Contract  No.  F44620- 
72-C-OO52,  31  December,  1975. 

Oberladstatter,  M.,  Ch.  Grissemann,  aind  Norbert  H.  Nessler. 

In  Situ  Measurements  of  the  Electric  Rock  Parameters 
by  Means  of  a Wenner  Four-Point  Arrangement , Scientific 
Report  No.  1,  AFCRL-TR-7^-0272 . Innsbruck,  Austria: 
Innsbruck  University,  Contract  No.  F44620-72-C-0052 , 

3 April  1974.  (AD  781  839) 

Porath,  H.,  "A  Review  of  the  Evidence  on  Low-Resistivity 

Layers  in  the  Earth's  Crust,"  Geophysical  Monograph 
Series . Volume  14.  The  Structure  and  Physical 
Properties  of  the  Earth  * s Crust,  edited  by  John  G. 
Heacock.  Washington,  D.C.:  American  Geophysical 
Union,  1971.  (AD  735  998) 

Rokityanskiy , I.,  "Deep  Electromagnetic  Soundings  in  the 
Ocean,"  Geomagnetism  and  Aeronomy , X( 4): 668-672 

(1970). 

Rotanova,  N.M.,  "Construction  of  a Chart  of  Conductivity  and 
of  the  Thickness  of  the  Nonconductive  Layer  from 
Geomagnetic  Data,"  Geomagnetism  and  Aeronomy , VI(1): 
89-92  ( 1966) . 

• "Distribution  of  the  Electric  Conductivity  of  the 

Earth,"  Geomagnetism  and  Aeronomy,  VII( 3) :435-438 

(1968). 

Spies,  Kenneth  P.  and  James  R.  Wait,  "Influence  of  the  Earth's 
Curvature  and  the  Terrestrial  Magnetic  Field  on  VLF 
Propagation,"  Journal  of  Geophysical  Research.  65: 

2325-2331  (196077^ 

Staras,  H.,  "Analysis  of  a Natural  VLF  Slot  Antenna,"  Symposium 
on  Electromagnetic  Theory  and  Antennas . Proceedings 
of  a Symposium  held  at  Copenhagen,  Denmark,  June 
1962,  edited  by  E.C.  Jordan.  New  York:  Pergamon 
Press,  Macmillan  Co.,  I963. 

Sutton,  G.H.,  G.L.  Maynard,  and  D.M.  Hussong,  "Widespread 
Occurance  of  a High-Velocity  Basal  Layer  in  the 
Pacific  Crust  Found  with  Repetitive  Sources  and 
Sonobuoys,"  Geophysical  Monograph  Series.  Volume 
14.  The  Structure  and  Physical  Properties  of  the 
Earth ' s Crust . edited  by  John  G.  Heacock.  Washing- 
ton, D.C.:  American  Geophysical  Union,  1971. 


68 


Troitskaya,  V.A.  and  O.M.  Barsukov, "Studies  in  Geomagnetism 
and  Aeronomy  State  of  Investigations  in  the  Fifeld  of 
Terrestrial  Currents,"  I ssledovaniva  Po  Geomagnet i zmy 
i Aeronomii.  Wright-Patterson  APB,  Ohio:  Foreign 
Technology  Division,  I963.  (FTD-TT-64-640)  (AD  610 

419) 

Wait,  James  R.,  "The  Theory  of  Electromagnetic  Surface  Wave 

of  Geologic  Conductivities,"  Geof isica  pura  e 

Annlicaccia.  (Milano),  28( 4?) :ri954) . 

, and  J.A.  Fuller.  Characteristics  of  an  Insulated 

Linear  Antenna  in  a Fluid  or  Air  Filled  Borehole. 
Technical  Report  No.  5*  Boulder,  Colorado:  Cooper- 
ative Institute  for  Research  in  Environmental 
Sciences,  University  of  Colorado  1 and  Institute  for 
Telecommunications  Sciences,  Office  of  Telecommuni- 
cations, Department  of  Commerce,  Contract  No.  NA- 
ONR-15-71,  28  February  1973.  (AD  757  977) 

• and  Kenneth  P.  Spies.  Magneto -Telluric  Fields  for 

a Segmented  Overburden . Technical  Report  No  7. 

Boulder,  Colorado:  Institute  for  Telecommunications 
Sciences,  Office  of  Telecommunications,  U.S.  Depart- 
ment of  Commerce,  Contract  No.  NA-ONR-I5-7I , 31  May 
1974.  (AD  781  740) 

, and  A.M.  Wahler.  On  the  Measurement  of  Ground 

Conductivity  at  VLF , NBS  Report  5037.  Boulder,  Colo- 
rado: Boulder  Laboratories,  National  Bureau  of  Stan- 
dards, U.S.  Department  of  Commerce,  14  December 
1956.  (AD  123  137) 

Wyllie,  P.J.,  "A  Discussion  of  Water  in  the  Crust,"  Geo- 
physical Monograph  Series , Volume  l4.  The  Structure 
and  Physical  Properties  of  the  Earth ' s Crust . edited 
by  John  G.  Heacock.  Washington  D.C.:  American 
Geophysical  Union,  1971. 


Electromagnetics 

Berdichevskiy,  M.,  L.L.  Van'yan,  and  E.B.  Faynberg,  "Theoret- 
ical Principles  in  Using  Electromagnetic  Variations 
to  Study  the  Electric  Conductivity  of  the  Earth," 
Geomagnetism  and  Aeronomy.  IX(  3) : 46 5-467. 

Berloffa,  E,  Dieter  Fritz,  R.  Klein,  aJ.  Theoretical  and 
Experimental  Studies  of  VLF  and  LF  Waves . Final 
Scientific  Report,  AFCRL  73-0174.  Innsbruck,  Aus- 
tria: Electronics  laboratory.  Physics  Institute, 
University  of  Innsbruck,  28  February  1973.  (AD  758 
778) 


69 


Bitterlich,  Wolfgang.  Measurement  of  the  Conductivity  and 

Dielectric  Constant  of  a Homogeneous  Medium  by  Means 
of  Wave  Propagation . Scientific  Report  No.  7»  AFCRL- 
ZE-OhkO . Innsbruck,  Austria:  VLF  Project,  25  July 
1968.  (AD  677  587) 

, On  the  Propagation  of  VLF  Waves  in  Solids , First 

Annual  Report,  AFCRL^2-378.  Contract  No.  61(052)- 
490,  December  I96I.  (AD  282  083) 

, On  the  Propagation  of  VLF  Waves  in  Solids , Second 

Technical  Annual  Summary  Report,  AFCRL-63-ioo . 
Innsbruck,  Austria;  Contract  No.  61(052) -490, 

20  December  1962.  (AD  402  3l4) 

, the  Propagation  of  VLF  Waves  in  Solids . Fourth 

Technical  Annual  Summary  Report,  AFCRL-65-203. 
Innsbruck,  Austria:  Contract  No.  61(052) -490, 

December  1964.  (AD  613  593) 

, On  the  Propagation  of  VLF  Waves  in  Solids . Final 

Scientific  Report,  AFCRL-66-211 . Innsbruck,  Austria: 
Contract  No.  61(052)-490,  December,  I965.  (AD  63I 
046). 

• On  the  Propagation  of  VLF  Waves  in  Solids  and  Liquids , 

Final  Scientific  Report,  AFCRL-69-OI52 . Innsbruck, 
Austria:  VLF  Project,  Contract  No.  61(052) -902, 

31  December  I968.  (AD  687  751) 

, VLF/LF  Antennas  and  Waves  in  Liquid  and  Solid  Media. 

Final  Scientific  Report,  AFCRL-7I-OI97 . Innsbruck, 
Austria:  VLF  Project,  Contract  No.  F6l( 052) -69-C- 
0007,  31  December  1970.  (AD  725  768) 

, and  Norbert  Nessler.  Propagation  of  Long  Electro- 
magnetic Waves  Through  Rock.  Scientific  Report  No.  2, 
AFCRL-67-O361 . Innsbruck,  Austria:  VLF  Project, 
Contract  No.  61(052) -902,  15  April  I967. 

, and  0,  Worz.  On  the  Frequency  Dependence  of  the 

Electric  Parameters  of  Rock.  Interim  Scientific 
Report  No.  15 1 AFCRL^5-686 . Innsbruck,  Austria: 
Contract  No.  6l(052)-490,  15  June  1965.  (AD  622  072) 

Bloomquist,  M.,  H.W.  Smith,  F.X.  Bostick,  ^.  , Propaga- 
tion Characteristics  of  Micropulsations . Final 
Report,  AFCRL-67-05^0 . Austin,  Texas:  Electrical 
Engineering  Research  Laboratory,  University  of 
Texas,  Contract  No.  19( 628) -4o46 , 1 August  1967. 

(AD  661  284) 


70 


1 


J 


Dolukhanov,  M.  "Underground  Propagation  of  Radio  Waves," 

Radio.  ( 1) :42-43( 1970) . In  Russian.  (JPRS  4989^) 

Grobner,  0.  and  W.  Bitterlich.  Propagation  of  VLF  Waves  in 
Highly  Conducting  Medium.  Technical  Note  7.  Inns- 
bruck, Austria:  Contract  No.  61(052) -490.  (AD  434 

776). 

Johler,  J.R.  and  R.L.  Lewis,  "Extra  Low  Frequency  Terrestrial 
Radio-Wave  Field  Calculations  with  the  Zonal  Har- 
monics Series,"  Journal  of  Geophysical  Research, 

74( 10) : 2459-2470  (15  May  I968) . (AD  69O  421) 

Kaspar,  M.  and  J.  Pecen,  "Detection  of  Caves  in  a Karst 
Formation  by  Means  of  Electromagnetic  Waves," 
Geophysical  Prospecting.  23( 4) ; 611-621 . 

Keller,  G.V.  Electrical  Properties  of  the  Earth ’ s Crust , 

Final  Scientific  Report,  for  I5  April-31  October 
1970*  Golden,  Colorado:  Colorado  School  of  Mines, 
Contract  No.  N00014-70-C-0290 , November  1971.  (AD 
739  739) 

, Recommendation  of  the  Working  Group  on 

Electromagnetic  Methods . Golden,  Colorado:  Colorado 
School  of  Mines,  1972.  (AD  740  122) 

Kellner,  W.  Electrical  Rock  Characteristics  From  the  Field 

Structure  of  a Magnetic  Dipole  Immersed  in  a Conduc- 
ting Medium.  Scientific  Report  No.  1,  AFCRL  -69-0281 . 
Innsbruck,  Austria:  VLF  Project,  Contract  No.  F 6I- 
O52-69-C-OOO7,  15  April  1969.  (AD  690  167). 

King,  R.W.P.  Antennas  In  Dissipative  Melia  and  in  Subsurface 
Communication . Final  Scientific  Report,  for  the 
period  1 July  1972-30  June  1975#  AFOSR-TR-1593. 
Cambridge,  Massachusetts:  Gordon  McKay  Laboratory, 
Harvard  University,  Contract  No.  F44620-72-C-002i , 

1 September  1975*  (AD  AOI9  682) 

Lukavec , R . Measured  Properties  of  Small  Self -Resonant 

Helical  Antennas  at  Low  Frequencies  for  Application 
in  Dissipative  Media.  Scientific  Report  No.  8,  AFCRL- 
69-0096 . Innsbruck,  Austria:  VLF  Project,  Contract 
No.  6l(052)-902,  1 December  I968.  (AD  684  578) 

l*«ssler,  Norbert  H.  Modelling  Studies . Technical  Note  1, 
Research  on  Propagation  of  VLF  Waves  in  Liquid  and 
Solid  Media.  Innsbruck,  Austria:  VLF  Project, 
Contract  No.  6l(052)-490,  1 July  I963.  (AD  418  336). 


71 


1 


A 


Nessler,  Norbert,  M.  Oberladstatter , and  Dieter  Fritz.  Fre- 
quency Dependence  of  Electrical  Rock  Parameters  in 
the  VLF  Range . Final  Report.  Innsbruck,  Austria: 
Electronics  Laboratory,  Innsbruck  University,  Contract 
No.  F44620-72-C-0052,  31  December  1975. 

Oberladstatter,  M.,  Ch.  Grissemann,  and  Norbert  Nessler.  In 
Situ  Measurements  of  the  Electric  Rock  Parameters  by 
Mpans  of  a Wenner~Four-PoTht  Arrangement . Scientific 
Report  No.  1,  AFCRL-TR-74-0272 . Innsbruck,  Austria: 
Innsbruck  University,  Contract  No.  F44620-72-C-0052 , 

3 April  1974.  (AD  781  839) 

O'Neal,  I.P.,  M.J.  McGowan,  F. Kelly,  and  F.J.  Rhoads,  VLF 
Propagation  Predictions  for  Disturbed  Conditions, 
Report  8017,  Naval  Research  Labs,  Electromagnetic 
Propagation  Division,  Communication  Sciences  Branch, 

27  July  1976. 

Ogorodneychuk,  I.F.,  I.  Ya.  Zhuravlev,  and  A. I.  Bykh,  "Remov- 
ing the  Electromagnetic  Field  of  a Low-Frequency 
Dipole  in  Rocks,"  Priborv  I^  Sistemy  Avtomatiki 
(Automatic  Devices  and  Systems) , No.  8:70-74,  1969. 

(AD  A008  265) 

Schlessinger , L.  and  W.J.  Karzas.  The  Fields  of  Buried  ELF 
Antennas . Topical  Report,  DNA  3545T . Santa  Monica, 
California;  R & D Associates,  Contract  No.  DNA  001- 
73-C-0237,  16  April  1975.  (AD  A012  685) 

Spies,  Kenneth  P.  and  James  R.  Wait.  Attenuation  Calculations 
for  a Two -Laver  Earth -Crust  Wave -Guide  for  Sub-Surface 
Electromagnetic  Propagation . Technical  Report  No.  2. 
Boulder,  Colorado:  Institute  for  Telecommunications 
Sciences,  Office  of  Telecommunications,  U.S.  Depart- 
ment of  Commerce,  Contract  No.  NA-ONR-15-71 . (AD  726 
75^) 

, "Influence  of  the  Earth's  Curvature  and  the  Terres- 
trial Magnetic  Field  on  VLF  Propagation,"  J oumal  of 
Geophysical  Research.  2325-2331  (i960). 

, Propagation  of  Electromagnetic  Waves  in  the  Earth 

Crust  Waveguide  From  ELF  to  VLF . Technical  Report  No. 
1,  Boulder  Colorado:  Institute  for  Telecommunication 
Sciences,  Office  of  Telecommunications,  U.S.  Depart- 
ment of  Commerce,  Contract  No.  NA-ONR-15-71*  1 May 
1972.  (AD  7U8  331) 

Tinhofer,  G.  Frequency  Dependence  of  the  Field  Strength 
Components  of  a Dipole  Embedded  in  a Conducting 
Medium.  Scientific  Report  No.  5*  AFCRL-68-0174. 


72 


Innsbruck,  Austria:  VLF  Project,  Contract  No.  61 
(052) -902,  30  November  I967.  (AD  668  518) 

Wait,  James  R.  "Analytical  Investigations  of  Electromagnetic 
Wave  Propagation  in  the  Earth's  Crust,"  Geophysical 
Monograph  Series , Volume  l4.  The  Structure  and 
Physical  Properties  of  the  Earth ’ s Crust , edited  by 
John  G.  Heacock.  Washington  D.C.;  American  Geophysi- 
cal Union,  1971. 

, "Electromagnetic  Propagation  in  an  Idealized  Earth 

Crust  Waveguide,"  Radio  Science . (New  Series),  l{S)i 
913-924(  August  1966T!  Cad  638  496) 

, Electromagnetic  Waves  in  Stratified  Media,  second 

edition.  New  York:  Pergamon  Press,  1970. 

• "Note  on  Calculations  of  Propagation  Parameters  for 

an  Idealized  Earth-Crust  Wave-Guide,"  Geophysical 
Monograph  Series . Volume  14,  The  Structure  and  Phys- 
ical Properties  of  the  Earth ' s Crust . edited  by  John 
G.  Heacock.  Washington,  D.C.:  American  Geophysical 
Union,  1971. 

, On  Mode  Conversion  in  the  Earth-Crust  Waveguide , 

Technical  Report  No,~4.  Boulder,  Colorado:  Cooper- 
ative Institute  for  Research  in  Environmental  Sciences, 
University  of  Colorado,  and  Institute  for  Telecommun- 
ications Sciences,  Office  of  Telecommunications,  U. 

S.  Department  of  Commerce,  Contract  No.  NA-ONR-15-71 1 
27  October  1972.  (AD  756  693) 

, "Propagation  of  Electromagnetic  Pulses  in  Terrestrial 

Waveguides,"  IEEE  Transactions  on  Antennas  ^d  Prop- 
agation. AP-13( 6) :901~918  (November  I965) . (AD  628 
240Ti 

, "Propagation  of  Pulses  in  Dispersive  Media,"  Scien- 
tific Report  No.  30,  AFCRL-66-435 1 Radio  Science 
Journal  of  Research.  NBS/USNC-URSI ."ggoTli) : 1387- 
l401  (November  19^5) . (AD  625” I39) 

. "Some  Factors  Concerning  Electromagnetic  Wave  Prop- 
agation in  the  Earth's  Crust,"  Proceedings  of  the 
IEEE.  ^(8)  1 1020  (1966). 

, "Terrestrial  Propagation  of  VLF  Radio  Waves  (A 

Theoretical  Investigation),"  Journal  of  Research 
of  the  National  Bureau  of  Standards . D.  Radio  Prop- 
agation. 64(  2)  ri63-203  (March -April  19601"^ 


73 


r 


Wait,  Jajnes  R.,  "The  Theory  of  Electromagnetic  Surface  Waves 
of  Geologic  Conductivities,"  Geof isica  pura  e Appli- 
ccia.  (Milano),  ^(47)j  ( 195^7"^ 

, Underground  Electromagnetic  Propagation , Preliminary 

Results  for  the  Period  1 August  1970-31  December 
1970,  ARPA  Order  No.  136I.  Boulder,  Colorado:  NOAA 
Environmental  Research  Laboratories;  and  Office  of 
Telecommunications,  Institute  for  Telecommunications 
Sciences,  Contract  No.  910  1468,  I5  January  1971. 

(AD  882  757) 

Viggh,  M.E.,  "Modes  in  Lossy  Stratified  Media  with  Application 
to  Underground  Propagation  of  Radio  Waves,"  IEEE 
Transactions  on  Antennas  and  Propagation . AP-ll(  3)  ’• 
318-323  (May  1963)7 


Communications 

Beitchman , J . G . , Soviet  Micropulsations  Research  for  Possible 
Communications  Application . Advanced  Projects  Research 
Agency,  ARPA  No.  2225,  21  February  1973.  (AD  759  380) 

i 

i Bitterlich,  Wolfgang.  On  the  Propagation  of  VLF  Waves  in 

I Solids.  First  Annual  Report,  AFCRL^2-378.  Contract 

I No.  61(052) -490,  December  I96I.  (AD  282  083) 

I 

1 , On  the  Propagation  of  VLF  Waves  in  Solide , Second 

I Technical  Annual  Summary  Report,  AFCRL-63-IOO . 

i Innsbruck,  Austria:  Contract  No.  61(052) -490, 

! 20  December  I962.  (AD  402  314) 

' Hansen,  R .C ., "Radiation  and  Reception  with  Buried  and  Sub- 

merged Antennas,"  IEEE  Transactions  on  Antennas  and 
Propagation.  AP-11(  Tm07-21 5 ( Mav  1961). 

, "Radiation  and  Reception  with  Buried  And  Submerged 

, Antennas,"  Electromagnetic  Theory  and  Antennas , 

edited  by  E.J.  Jordan.  Copenhagen,  Denmark:  Pro- 
, ceedings  of  a Symposium,  June  1962. 

i 

J Harwood,  J.  and  F.  Homer,  "Investigation  of  Atmospheric 

Radio  Noise  at  VLF,"  Proceedings  of  the  IEEE , 10 3B: 

743-751. 

King,  R.W.P.  Antennas  in  Dissipative  Media  and  in  Subsur- 
face Communic  ations . Final  Scientific  Report , for 
the  period  1 July  1972-30  June  1975 1 AFOSR-TR-75- 
1593.  Cambridge,  Massachusetts:  Gordon, McKay  Labor- 
atory, Harvard  University,  Contract  No.  F44620-72-C- 


74 


0021,  1 September  1975.  (AD  A019  682) 


Lukavec , R . Measured  Properties  of  Small  Self -Resonant  Hel- 
ical Antennas  at  Low  Frequencies  for  Application 
in  Dissipative  Media.  Scientific  Report  No.  8,  AFCRL- 
^-0096 . Innsbruck,  Austria:  VLF  Project,  Contract 
No.  61(052) -902,  1 December  1968.  (AD  684  578) 

Maxwell,  E.L.,  and  A.D.  Watt,  "Characteristics  of  Atmospheric 
Radio  Noise  from  1 to  100  KC , " Proceedings  of  the 
IRE.  ^*.787-794  (June  1957). 

Nessler,  Norbert  H.  and  W.  Bitterlich,  "Considerations  on 

Calibration  of  a Ferrite  Antenna  for  Measurements," 
Internationale  Elektronische  Rundschau . 25( 11) ; 277 
(1971)7 

Obayashi,  T.,  "Measured  Frequency  Spectra  of  VLF  Atmospherics," 

J oumal  of  Research  of  National  Bureau  of  Standards . 

D.  Radio  Propagation . 64d( 1) ; ^January -February  i960). 

Pierce,  E.T.  and  F.  Hepburn,  "Atmospherics  with  VLF  Components," 
Nature.  £2:837  (1953). 

Schlessinger , L.  and  W.J.  Karzas.  The  Fields  of  Buried  ELF 
Antennas . Topical  Report,  DNA  35^5T . Santa  Monica, 
California:  R & D Associates,  Contract  No.  DNA  001- 
73-C-0237,  16  April  1975.  (AD  A012  685) 

Staras,  H.  "Analysis  of  a Natural  VLF  Slot  Antenna,"  Sym- 
posium on  Electromagnetic  Theory  and  Antennas , 
Proceedings  of  a Symposium  held  at  Copenhagen,  Denmark, 
June  1962,  edited  by  E.C.  Jordan.  New  York:  Pergamon 
Press,  Macmillan  Co.,  I963. 

Tinhofer,  G.  Frequency  Dependence  of  the  Field  Strength 
Components  of  a Dipole  Embedded  in  a Conducting 
Medium,  Scientific  Report  No.  5,  AFCRL-68-0174. 
Innsbruck,  Austria:  VLF  Project;  Contract  No.  61(052)- 
902,  30  November  1967.  (AD  668  518) 

Wagner,  C.A.  and  M.  Balser,  "Measurement  of  the  Spectrum  of 
Radio  Noise  from  50  to  100  cycles/second,"  Journal 
of  Research  of  the  National  Bureau  of  Standards . D. 

Radio  Propagation.  6^4Dr 4)  :7uly-August  i960 . 

Wait,  James  R.  and  J.A.  Fuller.  Characteristics  of  an 

Insulated  Linear  Antenna  in  a Fluid  or  Air  Filled 
Borehole . Technical  Report  No.  5.  Boulder,  Colorado: 
Cooperative  Institute  for  Research  in  Environmental 
Sciences,  University  of  Colorado;  and  Institute  for 


75 


1 


Telecommunications  Sciences,  Office  of  Telecommun- 
ications, U.S.  Department  of  Commerce,  Contract  No. 
NA-ONR-I5-71 , 28  February  1973.  (AD  757  977) 

Wheeler,  H.A.,  "Useful  Radiation  From  an  Underground  Antenna," 
Journal  of  Research  of  the  National  Bureau  of  Stan- 
dards, 6 5D ( 1 ) : ( J anuary -F e br uary  19617"! 

Zhuralev , I.Ya.,  "Selection  of  Operating  Frequencies  for 

Underground  Telemetry,"  Priborv  i Sistemv  Avtomatiki. 
No.  5* 124-128( 1967) » in  Russian. 


76 


Appendix  As  List  of  Geologic  Terms 


Source*  Nelson,  A.  and  K.D.  Nelson.  Dictionary  of  Applied 
Geology.  Mining,  and  Cjyjl  Engineering.  New  York: 
Philosophical  Library,  Inc.,  I967. 


Silicate  of  calcium,  magnesium,  iron,  and  aluminum, 

(CaMgFeAl)2( AlSi)20g.  Composition  is  variable. 

Greenish  black  to  black  in  color. 

4 

A fine-grained,  often  partly  glassy,  basic  igneous  | 

rock  occuring  mainly  as  dykes,  sills,  and  lava  flows; 
essential  minerals  are  plagioclase  and  augite.  Dark 
grey  to  black  in  color. 

Basic  Rock.  Igneous  rock  containing  between  about  and 

45^  silica,  with  a high  proportion  of  dark  minerals 
and  no  free  quartz. 

Borehole . A hole  drilled  into  the  ground,  from  the  surface 
or  from  underground  workings,  to  secure  geological 
information;  also  used  for  the  extraction  of  water, 
oil,  natural  gas,  salt,  sulphur,  and  for  many  other 
mining  purposes.  A borehole  may  range  from  an  inch 
or  so  to  well  over  12  inches  in  diameter  and  drilled 
at  any  angle  to  many  thousands  of  feet  depth. 

Core.  The  cylindrical  piece  of  rock  obtained  during  dril- 
ling with  a diamond  drill  or  shot  drill.  The  core 
is  recovered  from  the  core  barrel  and  may  be  from 
less  than  an  inch  to  20  feet  or  so  in  length. 


Augite. 


Basalt. 
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Crystalline  Form.  A mineral  is  crystalline  when  the  atoms 


1 


1 
i 

are  arranged  in  a definite  geometrical  pattern.  | 

I 

There  are  six  principal  crystal  systems  and  these  j 

are  specified  by  the  relative  lengths  and  angular  j 

relationships  of  their  crystallographic  axes.  j 

Dunite . A variety  of  peridotite  consisting  of  almost  pure 
olivine . 

Dyke . A wall-like  mass  of  igneous  rock  that  fills  a fissure 
in  pre-existing  rocks;  often  form  a group  or  swami 
in  which  each  member  is  roughly  parallel. 

Feldspar.  A family  of  rock  forming  minerals  consisting  of 
aluminous  silicates  of  potassium,  sodium,  calcium, 
or  barium.  Greyish-white  to  red  in  color. 

Fissure . An  extensive  fracture,  crack,  or  break  in  rock 
formations. 

Format inn.  A group  of  rocks  which  have  certain  character- 
istics in  common  and  which  were  deposited  about  the 
same  geological  period  and  constitute  a convenient 
unit  for  description. 

Fracture . A break,  crack,  or  any  kind  of  discontinuity  in 

a mass  of  rock,  coal  or  orebody.  It  may  be  produced 
by  tensile  stress  or  shear  stress. 

Granite . An  important  and  widely  distributed  igneous  rock; 

coarse-grained,  mostly  deep-seated  or  plutonic.  It 
consists  of  quartz,  orthoclase  feldspar,  and  mica. 

Whitish-grey,  grey,  reddish  to  pink  in  color. 

78 

L 


1 


Igneous  Rock.  A rock  formed  by  the  solidification  of  molten 
matter  (magma)  from  the  earth's  interior. 

Lustre . The  appearance  of  the  surface  of  a mineral  in 

reflected  light.  The  lustre  may  differ  in  degree 
and  kind  on  different  crystal  surfaces. 

Magma.  Molten  masses  of  igneous  rock  before  they  have 

crystallized;  silica  is  the  dominant  component.  On 

cooling,  the  other  oxides  of  the  magma  combir  with  | 

the  silica  to  form  the  minerals  known  as  rock-forming  I 

silicates.  I 

i 


Metamorphic  Rock.  A sedimentary  or  igneous  rock  which  has 
been  altered  to  such  an  extent  as  to  erase  all 
resemblance  to  the  original  deposit. 

Mica.  Original  constituent  of  plutonic  and  volcanic  rocks; 
composed  of  complex  silicates  of  aluminum  and 
potassium.  Distinguished  by  its  splendent  pearly 
lustre,  softness,  and  brown-black  or  silvery -white 
color. 

Olivine.  Magnesium  iron  orthosilicate,  (M^e)2Si0j^,  in 

which  Fe  replaces  part  of  the  Mg.  Shades  of  green 
and  brown  in  color, 

Orthoclase.  Potassium  aluminum  silicate,  KAlSi^Og.  White 
to  pink  to  grey  in  color. 

Overburden.  The  barren  soil,  rock  or  waste  overlying  coal, 
or  an  orebody  at  a mine. 

Oxides.  Compounds  of  oxygen  with  another  element. 
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Peridotite.  A coarse-grained  ultrabasic  igneous  rock  composed 


dominantly  of  olivine.  Dark  colored. 

Plagioclase.  Silicates  of  aluminum  with  calcium  and/or 

sodium.  Colorless  to  white,  yellow,  red,  green, 
or  grey  in  color.  Important  mineral  in  many  plutonic 
and  volcanic  rocks. 

Plutonic  Rock.  A major  group  of  deep  seated,  coarse-grained 
igneous  rocks.  These  rocks  solidified  at  depth 
below  the  surface. 

Pyroxenes.  A group  of  silicate  minerals  which  similar  in 

certain  chermical  and  physical  characteristics.  The 
dark-colored  minerals  are  silicates  of  iron,  calcium, 
and  magnesium,  sometimes  with  aluminum. 

Quartz . Oxide  of  silicon,  Si02.  Colorless  to  white. 

Rock-forming  Minerals.  The  minerals  which  constitute  the 

treater  part  of  a rock  and  are  necessary  to  give  it 
its  character  as  a whole. 

Sedimentary  Rock.  A rock  produced  by  the  weathering  and 

denudation  of  pre-existing  rocks,  and  the  deposition 
of  the  sediment  in  water  or  on  land.  On  the  basis 
of  origin,  a sedimentary  rock  may  be  mechanically, 
chemically,  or  organically  formed. 

Seismic  Waye  Method.  A geophysical  prospecting  method  which 
measures  the  yelocities  of  shock  wayes  when  passing 
through  different  rocks.  The  survey  data  are 
interpreted  by  geophysicists  in  terms  of  the  nature 
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I and  structure  of  the  rock  formation  in  depth.  The 

i 


velocities  of  shock  waves  vary  with  the  elastic 
constants  and  densities  of  the  rocks  penetrated. 

The  method  is  used  to  locate  faults,  various  geo- 
logic structures,  and  assist  in  mining. 

Silica.  Silicon  dioxide,  Si02i  a colorless  or  white  mineral; 
hard  and  insoluble;  very  abundant  in  nature. 

Sill.  A sheet  of  igneous  rock  of  considerable  lateral 
extent  with  fairly  uniform  thickness. 

Stratification.  The  arrangement  of  the  deposited  sediments, 
and  eventually  the  rocks,  in  regular  layers  or  beds. 

Tectonics . The  physical  forces  which  are  operative  in  the 
earth's  crust  and  their  results;  closely  allied 
to  structural  geology. 

Ultrabasic  Rock.  An  igneous  rock  containing  less  than  about 
silica. 

Vitreous . A lustre  resembling  that  of  broken  glass.  Magma 
which  has  cooled  rather  rapidly  at  the  earth's 
surface  is  usually  vitreous  or  glassy,  or  partly 
glassy  and  partly  crystalline. 
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Appendix  B;  Computer  Methods 
Waveguide  Integration 

In  order  to  solve  Eq  (36)  for  the  phase  function  u, 
it  was  necessary  to  integrate  the  expression  on  the  right- 
hand  side.  A standard  library  integration  subroutine  was 
used,  while  the  main  program  and  the  supporting  subroutines 
aind  functions  were  original.  The  chart  on  the  next  page 
shows  the  general  relationships  between  the  various  program 
units,  although  it  is  not  intended  as  a formal  flow  chart. 
Program  listings  have  not  been  included,  as  they  were  impor- 
tant only  in  producing  the  data,  not  in  the  method  of  its 
production,  nor  in  determining  the  final  system  feasibility. 

Main  Program;  NIAM.  The  main  program  served  primarily 
to  input  initial  data  and  to  output  results.  Among  the 
values  input  by  NIAM  were  the  initial  parameters  for  the 
minimizer,  MINUM;  values  for  the  frequency,  intrinsic  dielec- 
tric constant,  magnetic  permeability,  and  the  depth  of  the 
waveguide;  the  initial  value  of  the  complex  propagation 
constant,  k,  and  estimated  upper  and  lower  limits  to  the 
possible  values  for  it;  and  algebraic  and  arithmetic  con- 
stants, calculated  from  physical  constants,  to  simplifiy 
later  calculations.  Many  of  these  parameters  and  constants 
were  passed  to  other  program  units  by  common  storage.  The 


main  program  then  called  MINUM.  When  MINUM  returned  with  the 
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best  values  of  the  k-parameters , NIAM  called  FUNP,  which 
calculated  the  best  value  of  k and  the  phase  profile  and 
printed  both  of  them  out 

Subroutine  MINUM.  This  library  routine  worked  in  terms 
of  two  parameters  that  varied  between  zero  and  one  as  the 
value  of  k varied  between  its  given  limits.  The  parameters 
were  set  to  0.5  in  NIAM,  in  order  to  start  k in  the  middle 
of  its  given  range.  FUN  was  then  called,  which  calculated 
the  variance  between  the  value  of  the  phase  function  at  the 
bottom  of  the  waveguide  and  the  value  given  by  the  boundary 
condition.  MINUM  attempted  to  minimize  this  variance  by 
changing  the  parameter  (and  thus  k)  values  in  random, 
gradient,  average,  and  jump  step  increments.  The  random 
steps  were  determined  by  a random  parameter  input  in  NIAM. 

The  last  random  numbers  used  were  output  for  use  in  the 
next  program  run,  in  order  to  assure  a new  set  of  truly 
random  numbers. 

Subroutine  FUN.  This  original  subroutine  performed 
several  functions.  First,  the  parameters  passed  from  MINUM 
were  converted  to  k-values  using 

K = KRL  + PM( 1)(KRU-KRL)  + i{KIL  + PM(2)(KIU  - KIL)}  (103) 

This  K-value  was  then  used  to  calculate  the  initial  value  of 
the  phase  function  at  the  upper  waveguide  boundary.  Then 
BLCKDQ,  the  integrator,  was  called,  which  returned  a final 
value  for  the  phase  function  at  the  lower  boundary.  FUN 
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recalculated  the  boundary  conditions  for  the  lower  boundary 
and  compared  this  value  of  the  pahse  function  with  that 
returned  by  BLCKDQj 

Im(YS  - US) ( 10^) 

Im(US)  + Re(US)  X 10“^/ 

This  fvinction  was  used,  rather  than  simply  subtracting  the 
real  and  imaginary  parts  of  YS  and  US,  because  of  the  differ- 
ence in  magnitude  between  the  real  and  imaginary  parts  of 
the  values.  It  was  necessary  to  normalize  the  parts  and 
insure  that  the  values  to  be  minimized  would  always  remain 
positive.  Thus  FUN  calculated  the  upper  and  lower  boundary 
conditions,  called  BLCKDQ , and  calculated  the  variance  of 
the  value  of  the  final  integrated  phase  function  from  the 
boundary  value. 

Subroutine  BLCDQ . This  library  subroutine  started 
with  the  initial  boundary  values,  calculated  in  FUN,  and 
solved  the  complex  first-order,  ordinary  differential  equa- 
tion down  to  the  lower  boundary  of  the  waveguide.  The 
integration  proceeded  from  the  initial  conditions  using  a 
starting  scheme  modelled  after  Picard's  method  of  successive 
substitutions  and  an  eighth-order  predictor-corrector 
formula  to  perform  the  integration.  Error  control  was 
effected  using  local  polynomial  fit  and  derivatives  deter- 
mined by  the  differential  equation.  Tolerances  were  relaxed 
as  the  solution  passed  a zero,  but  not  when  it  tended  towards 
one.  Thus  BLCKDQ  integrated  down  from  the  upper  boundary 
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condition  to  the  lower  one,  and  output  the  value  of  the 
phase  function  at  the  lower  boundary  condition. 

Function  FCN . BLCKDQ  used  this  original  function 
subroutine  to  define  the  function  to  be  integrated.  Because 
u,  the  phase  function,  was  complex,  the  problem  was  solved 
as  two  ordinary  differential  equations  corresponding  to  the 
real  and  imaginary  parts  of  u.  As  BLCKDQ  integrated  through 
the  waveguide , FCN  was  called  to  calculate  the  value  of 
du/dz  at  each  given  depth. 

Function  EP.  This  original  function  served  to  define 
the  complex  dielectric  constant,  €.  It  was  defined  separately 
so  that  the  various  conductivity  profiles  could  be  used  to 
obtain  different  permittivity  functions.  The  conductivity 
profiles  were  expressed  in  polynomial  form,  which  will  be 
discussed  below.  In  the  final  forms,  after  verification  of 
the  program  had  been  completed,  a logical  IF-statement  was 
added  to  define  the  discontinuity  at  the  upper  boundary.  '' 

Thus  for  the  expressions  defining  the  upper  boundary  conditions, 
a negative  argximent  was  used,  leading  to  a definition  of  the 
conductivity  in  terms  of  the  conductivity  of  seawater- saturated 
rock.  Then  for  any  zero  or  positive  depth,  the  conductivity 
function  was  defined  by  the  polynomial  form. 

Subroutine  LOOKAT.  Since  it  was  necessary  to  have  a 
mechanism  that  stopped  the  integration  by  BLCKDQ,  this 
routine  was  called  at  several  points  in  the  integration. 

LOOKAT  was  used  in  two  forms i a library  version  that  tested 
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whether  the  given  final  depth  value  had  been  reached,  and 
an  original  version  that  also  tested  whether  the  conductivity 
ahd  reached  a value  high  enough  to  prevent  wave  propagation 
within  the  guide  ( it  was  assumed  that  the  conductivity 
decreased  within  the  waveguide  region  before  becoming  large 
once  again  at  greater  depths) . The  latter  version  of  LOOKAT 
allowed  the  depth  of  the  waveguide  to  be  determined  mainly 
by  the  conductivity  profile. 

Subroutine  FUNP.  This  subroutine,  called  by  NIAM 
after  the  return  of  MINUM,  was  essentially  identical  to  FUN 
except  that  it  ran  only  once.  The  best  parameters  returned 
by  MINUM  were  converted  to  k- values  by  FUNP  and  printed. 

FUNP  also  called  BLCKDQ  in  a DO-loop,  so  that  the  integration 
was  performed  in  a series  of  sixty  equal  steps.  The  •values 
of  the  phase  function  and  depth  were  printed  at  each  step, 
providing  a phase  plot  of  the  impedance  function. 

Summary . A typical  program  cycle  started  with  IVINUM 
choosing  a set  of  parameters  for  k.  The  propagation  constant 
was  then  calculated  by  FUN,  and  BLCKDQ  was  called.  The 
fiinctjon  defined  by  FCN  was  integrated  using  the  initial 
conditions  from  FUN , in  accordance  with  the  conductivity- 
depth  profile  in  the  function  EP.  When  LOOKAT  determined 
that  the  bottom  of  the  waveguide  had  been  reached,  BLCKDQ 
retuined  the  final  phase  function  values  to  FUN.  The 
variance  from  lower  boundary  conditions  was  calculated  by 
FUN,  which  then  returned  to  MINUM  for  a new  set  of  k-parameters. 
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When  the  designated  number  of  trials  had  been  made , MTNUM 
called  FUNP,  which  calculated  and  printed  the  best  value  of 
k.  FUNP  also  printed  the  real  and  imaginary  parts  of  the 
phase  function  as  a function  of  depth.  By  adjusting  the 
upper  and  lower  limits  on  k,  in  successive  program  runs, 
the  variance  of  the  phase  function  from  the  boundary  con- 
ditions was  made  as  small  as  possible,  usually  on  the  order 
_2 

of  10  . This  "fine-tuning"  enabled  accurate  results  for 

the  best  value  of  k. 

Conductivity  Data;  Program  FIT 

The  conductivity  data  was  obtained  from  published 
curves  of  empirical  data,  in  the  form  of  a set  of  discrete 
values.  To  facilitate  input  to  the  function  EP,  it  was 
necessary  to  compute  a functional  form  to  fit  these  discrete 
values.  A fourth-order  polynomial 

Az^  + Bz^  + Cz^  + Dz  + E = ln<r  (105) 

was  used.  The  set  of  equations  in  this  form,  derived  from 
the  discrete  data  were  expressed  as  a matrix  equation 

7 • A = S ( 106) 

where  "K  was  the  coefficient  matrix,  which  was  unknown,  "Z 
was  a square,  symmetric,  positive-definite  matrix  of  the 
known  depth  values,  and  5 was  the  conductivity  matrix. 

A standard  library  subroutine,  MESCHP,  was  used  to  decompose 

88 


r 


7 by  the  method  of  Choleski  into  the  form 

^ ( 107) 

“T 

where  L was  the  transpose  matrix  of  U and  E was  lower-left 
triangular.  Thus 

• y = ^ ( 108) 

and 

• A = y ( 109) 

were  solved  by  forward  and  backward  substitution.  Thus  A 
was  obtained t which  contained  the  coefficients  used  to  form 
the  conducitivty-depth  function.  A set  of  plots  were  made 
for  each  set  of  data  and  the  corresponding  functions,  showing 
the  accuracy  of  the  functional  fit  to  the  actual  data.  These 
plots  are  in  the  Appendix  following. 
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Appendix  C>  Conductivity  Profiles 


Contained  herein  are  the  computer-generated  plots 
of  the  conductivity  data  taken  from  empirical  curves  in  the 
literature,  and  their  corresponding  curve-fit  functions. 
Except  for  the  first  two,  they  are  all  applicable  to  the 
crust  under  the  continents,  either  because  the  sources  stated 
that  they  were  measured  on  the  continental  crust,  or  because 
the  geophysical  discussion  accompanying  them  in  the  liter- 
ature indicated  a continental  crust.  These  are  included 
because  of  the  age  of  the  oceanic  profiles  available  and 
because  there  are  only  two  specifically  oceanic  profiles  that 
could  be  found. 

Figures  11  and  12 t Profiles  from  H.  Mott  and  A.W. 

Biggs  (Ref  3»  324-326).  These  are  composite  oceauiic  profiles, 

based  on  two  other  profiles,  obtained  by  different  methods; 

The  near-surface  values  for  both  were  taken  from  studies 

made  by  Garland  and  Webster  on  basement  rock  in  western 

Canada  (Ref  20t27).  The  method  used  involved  comparison 

of  magnetic  and  electric  field  components  of  the  natural 

/ 

geomagnetic  field.  From  10-60  km  the  data  was  based  on 
experimental  measurements  by  Hughes,  performed  on  olivines 
at  various  temperatures  (Ref  28).  The  curves  themselves 
required  an  assumed  linear  temperature -depth  gradient  of 
19°/km  (Ref  3«324). 
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Figure  13t  Profile  by  Keller  From  Gallawa  and  Haidle 
(Ref  Z’jtkk) . This  continental  profile  was  based  on  the  average 
of  well  log  measurements  at  two  locations,  separated  by  about 
fifty  miles  and  in  differenct  rock  formations.  The  lower 
part  was  derived  from  in-laboratory  dry  rock  sample  measure- 
ments. A thermal  gradient  of  l6°/km  in  the  crust  was  assumed. 
As  this  profile  was  originally  contained  in  a private  commun- 
ication, it  was  first  published  by  Gallawa  and  Haidle,  even 
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Figures  and  15 » Profiles  by  Simmons,  Housely,  and 
Forin  (Ref  29)  from  Field  and  Pore  (Ref  4; 10-12).  These  two 
continental  profiles  were  the  results  of  laboratory  measure- 
ments based  on  two  relatively  new  and  controversial  concepts: 
1)  that  fluid-filled  microcracks  are  present  in  rock  removed 
from  deep  in  the  earth,  and  cannot  be  closed  by  reapplying 
pressures  and  temperatures  comparable  to  those  in  the  rocks’ 
original  locations.  Thus  laboratory  measurements  resulted 
in  significantly  higher  conductivities  than  may  be  present 
within  the  earth.  2)  The  partial  pressure  of  oxygen, 
different  for  deep  in  the  earth  from  at  the  surface,  has 
an  effect  on  laboratory  conductivity  measuremen,ts.  The  two 
profiles  vary  due  to  different  thermal  gradients:  BR, 

2k°/km;  BUS,  l4°/km. 
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Figure  l6t  Profile  from  Schwering.  Peterson,  and 
Levin  (Ref  30) ■ in  Levin  (Ref  21; 334).  This  continental 
profile  represents  the  concept  of  a two-part  conductivity 
function.  The  upper  segment  has  a decreasing  conductivity 
as  the  water  content  decreases  in  the  rock,  due  to  tempera- 
ture and  pressure  effects.  This  is  the  reason  for  referring 
to  the  profile  as  "dry" . The  lower  segment , below  ten 
kilometers,  reflects  the  increasing  temperature  effects, 
which  produce  dominant  electronic  conductivity  effects. 
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Figure  17 » Profile  from  Brace  (Ref  31»10)  in  Levin 
(Ref  21 t 334) . This  continental  profile  was  based  on  the 
concept  that  the  electrolytic  conductivity  from  interstitial 
fluids  was  the  primary  effect  down  to  a depth  of  30  km, 
where  the  electronic  effects  became  dominant.  The  presence 
of  a water  content  throughout  the  upper  segment  of  the  pro- 
file accounts  for  its  being  called  "wet". 
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The  assumption  was  made  that  the  attenuation,  as  a 
function  of  the  logarithm  of  the  frequency,  could  be 
expressed  as  a quadratic  equation: 


(P  - Pq)^  = ^a(a  - a^) 

(110) 

where 

P = lOgj^QW 

(111) 

and  where  a is  a constant  and  the  minimum  point  of 

the 

semilog  parabola  is  at  (Pgf  method  used 

late  the  minimum  frequency  was  as  follows. 

to  calcu- 

Given 

(p2,a2)  (P^,a^) 

(112) 

calculate 

~ '*■  ^3^ 

(113) 

Zj  ' *2(^1  + P3) 

(114) 

Zj  = ttjOj  + Bj) 

(115) 

1 

Th6  minimum  is 


(117) 


and 


e -B.  ^ ■ 


CSi32], 


(118) 


The  miniiniiin  attenuation  was  obtained  from  the  parabolic 
equation,  Eq  (110),  above.  The  minimum  value  was 
substituted  from  the  first  calculation,  along  with  two 
sets  of  values  from  Table  II . The  constant  a was  then 
eliminated  from  the  two  equations,  and  the  minimum  value 
was  calculated. 
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